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floating  probe  the  ion  content  flowing  from  a  microwave  discharge  of  hydrogen 
was  found  to  be  of  the  total  flow,  thus  ion  molecule  reactions  are 

unimportant^  Using  pressure  profiles  for  the  pulsed  reactant  it  was  calcula¬ 
ted  that  the  residence  time  in  the  reactor  is  approximately  50  milleseconds; 
for  CM  experiments  radiational  relaxation  of  HF  is  important  and  for  arrested 
relaxation  (steady  state  experiments)  a  background  of  relaxed  HF  may  be  a 
source  for  collisional  relaxation. 

The  interpretation  of  CM  experiments  requires  realistic  modeling  to  decon- 
volute  the  experimental  data.  Primary  processes  include:  initial  formation 
of  products,  collisional  and/or  radiative  relaxation,  and  physical  removal  of 
reactants  and  products  from  the  observation  volume.  Using  the  well  studied 
reaction,  F  +  H2  -->  HF ( v)  +  H,  as  a  bencnmark,  experiments  and  computer 
simulations  reveal  two  types  of  CM  experiments:  CM-PP  in  which  there  is  a 
repetitive  pressure  pulse  of  the  diatomic  species  and  CM-LP  in  which  the 
atomic  species  is  formed  repetitively  in  a  constant  pressure  environment. 

The  results  of  the  CM-PP  computer  simulations  are  understood  in  terms  of 
deviations  from  a  steady  state  system  and  by  the  relative  fluxes  of  the 
primary  processes  as  a  function  of  time.  CM-LP  experiments  are  easier  to 
analyze  than  CM-PP  experiments  because  of  pseudo  first  order  conditions. 

'Both  types  of  CM  experiments  provide  nascent  vibrational  energy  distributions. 
In  CM-PP  experiments,  obtaining  these  distributions  is  simplified  for  condi¬ 
tions  of  low  flow  (F  »  H2),  observation  times  must  be  short  enough  to  resolve 
the  fastest  relaxation  process,  and  the  reaction  vessel  must  be  evacuated 
between  pulse  sequences.  The  previous  CM-PP  experiments  in  which  the  nascent 
distribution  for  the  D  +  F2  was  reported  did  not  meet  these  requ i rements ;  the 
pulse  time  was  reduced  from  10  to  0.5  msec. 

CM-LP  experiments  also  provide  total  reaction  rate  constants  and  microscopic 
relaxation  rate  constants.  To  process  the  experimental  data,  an  iterative 
method  of  extracting  reaction  and  relaxation  rate  constants  for  vibrational 
levels,  (v),  below  the  maximum  vibrational  level  populated,  was  developed. 
Input  to  the  v-2  level  from  relaxation  of  the  v  level  was  decoupled  from  the 
reaction  input  to  the  v-1  level.  It  is  shown  that  reaction  and  deactivation 
rate  constants  obtained  by  analysis  of  specific  vibrational  levels  are  more 
accurate  than  those  obtained  by  analyzing  total  intensity-time  profiles. 

In  addition  to  the  F  +  benchmark  reactions,  experimental  results  for 

H(D)  +  F2,  F  +  HBr,  and  F  +  CH2C12  are  discussed  and  compared  to  other 
studies  when  available..  Discrepancies  between  the  previous  CM  results  and 
those  reported  here  can\be  resolved  by  the  high  reagent  flows  and  long 
observation  times  used  in  the  earlier  experiments. 
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I.  INTRODUCTION 


The  D2/F2  chemical  laser  is  comprised  of  elementary  chemical 
reactions  which  include  various  feedback  loops1-3  so  that  a  chain 
reaction  is  possible.  The  chain  is  propagated  with  the  production  of 
either  deuterium  or  fluorine  atoms.  For  typical  DF  lasers  the  chain 
length  is  estimated  to  be  about  30.  The  overall  laser  efficiency  is 
directly  related  to  the  production  of  DF  from  the  elementary  processes; 
an  increase  in  the  concentration  of  D  or  F  will  directly  enhance  the  DF 
production. 

A  simplified  laser  model  consists  of  DF  production  steps,  DF 
relaxation  steps  and  the  stimulated  emission  of  photons  (lasing).  This 
model  is  depicted  by 


The  laser  output  energy  can  be  optimized  by  increasing  the  production 
rates  and  decreasing  the  relaxation  rates.  In  order  for  this  to  be 
accomplished  detailed  knowledge  of  the  important  production  and 
relaxation  steps  involving  DF  must  be  known. 

To  display  the  complex  set  of  reactions  with  the  various  feedback 
loops  a  schematic  representation  is  shown  in  Figure  1-1.  The 
terminology  for  this  diagram  is  that  the  initial  state  is  defined  at  the 
tail  end  of  the  arrow  while  the  final  state  for  a  particular  process  is 
at  the  head  of  the  arrow.  This  representation  is  analogous  to  that  used 
in  models  demonstrating  transport  in  biochemical  systems.  Although  the 
diagram  is  complex,  it  is  easily  seen  that  any  process  that  increases 
DF  will  enhance  the  probability  of  lasing.  The  traditional 
representation  o£  these  processes  is  listed  in  Table  1-1.  The  net 
production  of  DF  can  be  placed  into  three  elementary  categories: 

1.  initial  (nascent)  production  (processes  IF  and  ID) 

2.  relaxation 

a.  collisional  (processes  2D,  2F,  2D2,  2F2  and  2M) 

a.  radiative  (processes  2  and  2hv) 

3.  enhancements  ± 

a.  reactant  translation  (processes  and  l£*) 

b.  reactant  vibration  (processes  IF  and  ID  ) 

c.  chain  branching  (processes  3DF  and  3D2) 

The  radiative  relaxation  processes  (2  and  2hv),  corresponding  to 
spontaneous  and  stimulated  emission,  are  well  understood?  their 
contributions  are  calculated  from  the  well  known  Einstein  coefficients4. 
The  enhancement  reactions  are  based  on  increased  reactivity  with 
reactant  energy  (3a  and  3b)  and  an  increase  in  chain  carrier  (3c) 
concentrations.  These  lead  to  thermal  and  population  explosions, 
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Table  1-1 

Elementary  Processes  Involved  in  the  D2/F2  Chemical  Laser 


Process 


Label+ 


PRODUCTION  (1) 


RELAXATION  (2) 
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ENHANCEMENTS 
EXCITED  REACTANTS 


F*  +  Do  --- >  DF*  +  D 


D*  +  F2  — >  DF  +  F 


IF* 

ID* 


F 

D 


+  D  * - >  DF*  +  D 

+  F2* — >  DF*  +  F 


IF** 

ID** 


CHAIN  BRANCHING  (3) 


DF*  +  F2  — ->  DF*  +  F  +  F 


D2*  +  F2  — ->  DF*  +  F  +  D 


3DF 

3D0 


— >  d2  +  f  +  f 


+  The  reaction  labels  are  generated  by  the  type  of  process 
(1,2,3)  followed  by  the  atomic  species  for  reaction  or  the 
deactivator  for  relaxation.  A  single  *  represents  translational 
excitation  while  a  **  corresponds  to  vibrational  excitation  in 
the  diatomic  reactant. 
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respectively 


Although  the  above  description  gives  an  overall  model  for  the  gross 
features  of  the  D5/F2  chemical  laser  system  a  more  detailed  model 
including  particular  vibrational  and  rotational  states  of  DF  (v  &  J, 
respectively)  must  be  invoked  to  understand  specifics  such  as  laser 
emission  wavelength/  output  power  etc.  Hart5  has  reported  on  the 
computer  model  for  the  DF  chemical  laser  and  the  need  for  an  enlarged 
and  reliable  data  base  for  the  rate  constants  of  the  elementary 
processes. 

The  partitioning  of  the  energy  released  in  the  DF  production 
reactions  (IF  and  ID)  have  ±>een  experimentally  characterized  by 
replacing  D-,  and  D  with  H^6-8  and  H9-1  ,  respectively.  Using  similar 
techniques  the  energy  distributions  for  DF  resulting  from  reactions  IF6- 
9,14  and  ID13  14  have  been  reported.  All  of  these  studies  were 
performed  with  the  reactants  in  a  vessel  at  room  temperature  or  below. 
Although  these  studies  provide  fundamental  information  for  the 
understanding  of  chemical  dynamics  they  may  not  be  appropriate  for 
modeling  the  D2/F2  chemical  laser.  High  temperatures  will  be  achieved 
since  reactions  IF  and  ID  are  highly  exoergic  and  the  D2/F2  reactor  is 
thermally  isolated.  These  high  temperatures  will  enhance  the  reaction 
rate  constants  and  effect  the  vibrational  energy  distribution  of  the  DF 
product.  Thus  in  order  to  understand  the  overall  operation  of  the  D2/F2 
chemical  laser  the  reaction  dynamics  of  the  fundamental  processes  must 
be  understood,  an  understanding  that  is  necessary  for  optimizing  the 
efficiency  of  these  chemical  lasers. 

Various  experimental  techniques  have  been  used  to  obtain  the 
nascent  vibrational-rotational  energy  distribution.  Host  of  these 
techniques  use  a  microwave  discharge  to  produce  the  atomic  species. 
Wang15  has  suggested  that  both  atoms  and  ions  (H  ,  D~,  D2-,  H2-,  etc) 
are  produced  in  the  atomic  source  and  the  possiblity  of  negative  ion- 
molecule  reactions  (NIMR)  must  be  considered  when  analyzing  earlier 
results.  Similar  NIMR's  have  been  observed  by  Leone  and  coworkers36. 
Wang15  estimates  the  rate  constants  for  NIMR's  to  be  100-1000  times 
faster  than  the  corresponding  atom  molecule  reaction.  Thus  an  ion 
"impurity"  (contaminant)  of  0.1%  would  allow  NIMR's  to  compete  favorably 
with  the  expected  neutral  reactions.  Measurements  pertaining  to  the  ion 
content  of  gases  flowing  thru  a  microwave  discharge  cavity  are  presented 
in  section  II  of  this  report.  These  results  indicate  that  NIMR's  are 
not  important  for  the  experimental  conditions  used  in  producing  D,  F  or 
H  atoms. 

Measurements  of  the  formation  and  relaxation  rate  constants  and  of 
nascent  distributions  have  been  obtained  using  steady  state  experiments 
and  monitoring  the  spectrally  resolved  infrared  emission 
(chemiluminescence).  These  steady  state  experiments  are  somewhat 
limited  by  the  relative  magnitudes  of  the  rate  constants  for  relaxation, 
reaction  and  infrared  emission  and  the  large  pumping  speed  necessary  to 
maintain  a  short  lifetime  for  the  excited  molecule.  Chemiluminescent 
Mapping  (CM)14'1  ,  a  pulsed  technique,  has  been  developed  as  an 
alternative  method  in  obtaining  formation  and  relaxation  information  for 
exoergic  bimolecular  exchange  reactions.  Both  time  and  spectral 
resolution  are  obtained  simultaneously. 


Initial  CM  results  for  reactions  IF  and  ID  have  been  reported. 
Although  the  technique  was  benchmarked  with  the  F  atom  reactions  the  F2 
results  were  suspicious.  In  addition  to  the  possibility  of  NIMR's  as 
discussed  above  the  primary  concern  was  with  the  'unknowns'  associated 
with  a  new  experimental  technique.  Section  II  of  this  report  includes 
diagnostics  of  the  CM  apparatus  while  section  III  presents  model 
calculations  which  simulate  the  CM  experiments.  These  simulations 
provide  the  means  for  analysis  of  the  experiments  as  well  as  providing 
the  optimum  experimental  conditions  for  abstracting  the  fundamental 
information  from  the  raw  experimental  data.  The  CM  experimental 
apparatus  and  data  processing  are  described  in  section  IV  while  the 
initial  experimental  results  from  CM  experiments  are  presented  in 
section  V.  These  results  provide  the  foundation  for  future  studies 
related  to  the  processes  involved  in  the  D2/ F2  chemical  laser. 
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II.  DIAGNOSTICS  OF  CM  APPARATUS 


Earlier1^  we  reported  nascent  vibrational  energy  distributions  for 
the  reactions 

H(D)  +  F 2 - >  HFV  (DFV)  +  F  vmax=<10{15)  II-l) 

F  +  H2  (D  2 ) - >HFV  (DFV)  +  H  (D)  vmax<=3  (4)  II-2) 

using  Chemiluminescence  Mapping  (CM)  in  which  the  diatomic  reactant  is 
pulsed  into  a  low  pressure  of  atomic  reactants.  A  specially  designed 
Evenson  type  microwave  cavity  was  used  to  generate  either  F,  D  or  H 
atoms.  The  CM  product  energy  distribution  for  the  F  atom  reactions  were 
in  agreement  with  those  reported  by  other  techniques,  however  the  F2 
results  indicated  a  higher  level  of  vibrational  excitation  than  that 
which  was  observed  in  steady  state  experiments. 

Although  the  technique  was  benchmarked  with  the  F  atom  reactions 
the  F2  results  were  suspicious.  Of  primary  concern  was  the  possiblity 
of  negative  ion-molecule  reactions  (NIMR);  similar  NIMR's  have  been 
observed  by  Leone  and  coworkers16.  It  has  been  suggested  by  Wang15  that 
earlier  chemiluminescence  results  in  which  a  microwave  discharge  was 
used  for  producing  the  atomic  species  are  due  to  NIMR's  and  not  from  the 
assumed  neutral  reactions.  The  postulated  negative  ions  (H“,  D~,  D2“, 
H?-,  etc)  are  allegedly  generated  by  the  microwave  discharge  source 
which  is  used  to  produce  atomic  species  (H,  D,  F).  The  negative  ion 
molecule  reactions  suggested  are: 

d-  +  F2 >DFv  +  F  +  e  »„,=<!!  1 1-3 ) 

F  +  d2-  — >  DFV  +  D-  vmax  =  <1^  H'4’ 

Wang15  estimates  the  rate  constants  for  reaction  II-3  to  be  1000  times 
faster  than  II-2  and  reaction  II-4  to  be  100  times  faster  than  II-2. 

Thus  an  ion  "impurity"  (contaminant)  of  0.1%  would  allow  reactions  II-3 
and  II-4  to  compete  favorably  with  II-l  and  II-2.  For  normal 
experimental  conditions,  the  reaction  involving  D2“  will  not  be 
important  since  its  lifetime  is  less  than  a  microsecond1'. 

Diagnostic  experiments  (unpublished)  in  our  laboratory  using  a 
single  probe  indicated  the  ion  concentration  from  a  microwave  discharge 
cavity  to  be  <  0.0003  (0.03  %)  of  the  atom  concentration.  If  the  NIMR's 
have  a  reaction  cross  section  100-1000  times  greater  than  the  neutral 
reaction  then  the  NIMR  products  could  account  for  up  to  30%  of  the 
observed  HF  (DF).  Since  there  are  inherent  problems  with  the  single 
probe  technique  it  appeared  conceivable  that  NIMR  may  be  interfering 
with  the  reactions  which  had  bee.ri  reported.  Preliminary  double  floating 
probe  experiments  were  reported16  with  a  comparable  upper  ion 
concentration  limit.  Hence  quantitative  analyses  of  ions  in  the 
discharge  region  and  in  the  reaction  vessel  with  a  double  floating  probe 
were  needed.  Additionally  the  pressure-time  profile  of  the  pulsed 
diatomic  molecule  is  necessary  in  order  to  adequately  model  the  CM 
results  and  abstract  the  nascent  population  distribution.  This  section 
includes  the  diagnostics  of  the  CM  apparatus;  the  results  can  be  applied 
to  any  system  in  which  atoms  are  formed  by  a  microwave  discharge  and/or 


a  similarly  pumped  reactor  (as  is  used  in  the  Arrested  Relaxation 
experiments:  see  reference  14  for  additional  citations) 


ION  CONTENT:  Double  Floating  Probe  Experiments 
Background: 

The  microwave  induced  plasma  (MIP)  which  is  formed  in  the  discharge 
region  of  the  microwave  cavity  has  been  reviewed19;  aside  from  atom 
production  the  MIP  hAf  been  used  successfully  as  a  source  for 
spectroscopic  studies20.  Although  the  overall  mechanism  involved  in  a 
plasma  is  not  completely  understood19  it  is  believed  to  primarily 
involve  the  ionization  of  the  gas  molecules  by  high  energy  electrons; 
metastables  producing  Penning  ionization  may  also  be  involved. 

The  use  of  a  double  floating  probe  technique  in  the  CM  system  was 
needed  to  reduce  the  disturbance  (withdrawal  of  charge)  of  the  plasma  by 
the  probe  itself  (especially  in  the  reaction  vessel  measurements)  and  to 
provide  a  reference  potential  (a  defined  space  potential)  which  is  not 
present  in  the  single  probe  technique.  The  double  floating  probe 
technique  has  been  applied  successfully  to  MIP's  by  Busch  and  Vickers21 
and  by  Brassem  et  al22  25 ;  the  characteristic  curve  (current  vs  applied 
voltage)  is  used  to  calculate  the  electron  concentration  and 
temperature.  Ideally  the  characteristic  should  be  symmetric  about  the 
origin  and  exhibit  saturation  currents.  The  non-ideal  characteristic  is 
due  to  probe  asymmetry  (different  surface  areas  or  collecting 
efficiencies),  different  floating  potentials  resulting  from  a  potential 
gradient  in  the  plasma  or  a  difference  in  contact  potential  between  each 
probe  and  the  plasma  (this  could  be  due  to  contamination). 

Experimental  Setup: 

The  specifics  of  the  CM  apparatus  have  been  reported  earlier14,26. 
The  microwave  cavity  encapsulated  a  12  mm  OD  quartz  flow  tube  leading 
into  a  continually  pumped  10  1  reaction  vessel  (<0.0001  torr);  a  1  mm 
restriction  at  the  end  of  the  tube  maintained  a  steady  pressure  of  1 
torr  in  the  discharge  region.  A  Raytheon  Model  PGM  10X2  microwave 
generator  with  variable  output  power  up  to  120  watts  was  used. 

In  order  to  assure  that  systematic  errors  were  absent  in  these 
results  a  set  of  benchmark  experiments  were  performed.  These 
experiments  probed  the  discharge  region  of  the.  MlP^and  can  be  directly 
compared  to  results  from  similar  MIP  systems21'23-25.  After  confirming 
the  benchmark  experiments  the  reaction  vessel  was  probed  for  ion 
content. 

The  experimental  setup  for  the  benchmark  experiments  is  shown  in 
Figure  II-l.  The  probes  were  designed  according  to  the  limitations  on 
probe  radius  and  probe  separation  assumed  in  the  theory  as  delineated  by 
Schott27.  Two  sets  of  cylindrical  probes  were  used:  copper  (0.10  +0.05 
mm  OD  by  2.00  +  0.05  mm  in  length)  and  tungsten  (1.00  +0.05  mm  OD  by 
2.00  +0.05  mm  in  length).  The  copper  probe  has  a  total  surface  area  of 
6.2  *  10-3  cm2  while  that  of  the  tungsten  probe  was  7.4  *  10“ 2  cm2;  the 


Reaction  Vessel 


Figure  1 1  - 1 .  Schematic  Diagram  for  Setup  of  Discharge  Tube  Experiments. 


*.  \ 


inactive  surface  area  was  insulated  with  "Corona  Dope"  (GC  Electronics), 
a  high  voltage  insulating  paint. 

The  probes  were  mounted  on  a  1/4"  plate  of  plexiglass  using  high 
vacuum  epoxy;  the  separation  was  6.0  and  2.90  mm  for  the  copper  and 
tungsten  probes,  respectively.  This  plate  was  attached  to  the  discharge 
tube  with  a  vacuum  seal  comprised  of  a  #9  glass  joint  and  'O'  ring.  The 
copper  probes  penetrated  2  mm  into  the  discharge  tube  while  the  tungsten 
probe  was  "1.4  mm  from  the  edge  of  the  discharge  tube.  Due  to  the 
plexiglass  mounting,  a  high  voltage  degassing  heater  could  not  be  used; 
unfortunately  this  enhances  the  chance  of  probe  contamination. 
Nonetheless  contamination  was  only  observed  near  the  probe  where  carrier 
concentrations  are  highest.  The  probe  was  nominal lv  degassed  with  a 
Tesla  coil  prior  to  measurments  as  reported  earlier28. 

An  Evenson  microwave  cavity2^  with  a  removable  end  cap  allowed  for 
easy  movement  so  that  the  distance  between  probes  and  microwave  cavity 
could  be  systematically  varied;  this  distance  (D)  was  measured  from  the 
probe  arms  to  the  bottom  of  the  microwave  cavity;  the  distance  from 
probe  to  the  center  of  cavity  is  0.75"  larger.  Alternatively,  one  end 
of  the  quartz  discharge  tube  was  connected  to  a  conventional  low 
pressure  gas  handling  system  and  the  other  end  to  a  3"  oil  diffusion 
pump. 

The  probes  were  connected  to  a  variable  voltage  power  supply;  a 
Keithly  417K  picoammeter  was  connected  in  series  with  one  of  the  probes, 
while  a  Dana  digital  voltmeter  was  connected  in  parallel.  All  of  these 
devices  were  isolated  from  external  ground  as  is  required  for  the  double 
floating  probe  technique. 

The  probes  for  the  reaction  vessel  experiments  were  made  from  the 
inner  elements  of  a  vacuum  tube.  The  tungsten  probe  had  an  inner  total 
surface  area  of  1.24  cm2  and  a  total  outer  area  of  34.3  cm2;  probe 
separation  was  1.78  cm.  Since  this  asymmetric  probe  has  a  surface  area 
ratio  <300  it  fulfills  the  requirements  for  a  floating  probe30.  This 
probe  was  mounted  on  a  glass  rod  which  in  turn  was  mounted  on  a  vacuum 
feed  through;  this  allowed  the  probe  to  be  adjusted  externally  without 
breaking  vacuum. 

Figure  II-2  is  the  characteristic  for  the  tungsten  probes  1.5"  (D) 
from  the  cavity;  this  curve  has  high  symmetry  and  passes  through  the 
origin.  The  large  sloping  saturation  regions  are  typical  of 
measurements  made  inside  the  discharge  tube. 

Characteristic  curves  were  determined  for  the  copper  probes  1.5" 
from  the  cavity  for  3  different  levels  of  microwave  power.  Although  the 
distance  from  the  cavity  is  the  same  as  that  for  the  tungsten  probes 
these  characteristics  show  an  increasing  deviation  from  the  origin  as 
the  power  is  increased  because  the  potential  gradient  across  the  probe 
increases  with  power.  The  large  saturation  slopes  also  indicate  the 
influence  of  the  positive  ion  species  at  high  voltage. 

Characteristics  obtained  for  the  probes  located  in  the  reaction 
vessel  at  60  and  96%  power  exhibit  large  Vapp  for  1^  =  0  ^ue  t0  probe 
asymmetry.  As  observed  for  the  discharge  tux>e  experiments  the  ion 
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Figure  II-2.  Double  floating  probe  characteristics:  data  taken  from 
experiment  #3-1  with  tungsten  probes  1.5"  from  cavity  and  60  watts 
microwave  power.  Ipj,  Ip  and  I  _  are  the  currents  at  probes  1  and 
2  due  to  ions  and  electrons,  respectively. 
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current  also  increases  with  microwave  power. 

The  characteristic  for  the  tungsten  probes  near  the  discharge 
region  did  not  exhibit  plateaus;  this  was  a  diagnostic  indicating  the 
presence  of  probe  contamination.  Character stic  curves  3  ike  these  were 
not  used  in  calculating  ion  densities.  The  shape  of  these  curves  are 
similar  to  those  of  Thornton31. 

RESULTS; 

Electron  temperature  can  be  determined  by  the  logarithmic  plot, 
equivalent  resistance  and  intercept  methods;  the  difference  in 
calculated  temperature  between  these  methods  is  <5%  2  .  The  equivalent 
resistance  and  intercept  methods  were  used  in  this  work.  The  electron 
temperature  does  not  depend  on  the  probe  area,  contact  potential  or 
difference  in  plasma  potential  between  probes.  The  electron 
concentration  was  calculated  from  the  equation  derived  by  Massey32: 

isat=0  *  5*N0*A*e* (kTe/m) 1/2 

where  N0  is  the  number  of  positive  ions/cc,  A  is  the  surface  area  of  the 
probes,  Te  is  the  electron  temperature  and  e  is  the  electronic  charge. 
This  equation  has  been  derived  for  spherical  probes;  for  cylindrical 
probes  only  an  order  of  magnitude  concentration  can  be  extracted.  The 
results  for  the  double  floating  probe  experiments  (complete 
characteristics)  are  tabulated  in  Table  II-l;  the  discharge  tube  and 
reaction  vessel  results  will  be  discussed  separately. 


Discharge  Tube  Results 

2  S  12 

Brassem  has  measured  electron  densities  of  approximately  2  *  lO-1^ 

per  cc  and  electron  temperatures  of  ”  100,000  K  under  similar  conditions 
to  these  experiments  inside  the  microwave  cavity.  These  values  are 
larger  than  experiments  #7-1  and  8-1  of  the  present  work  which  are  at 
lower  microwave  power  and  not  in  the  microwave  cavity  (i.e.  0.75" 
removed).  As  seen  in  experiments  #1-1  and  4-1  or  2-1  and  5-1  or  8-1, 9-1 
and  11-1  the  ion  densities  decrease  dramatically  with  distance.  For  an 
atmospheric  argon  plasma  Goode  and  Otto33  showed  a  decrease  of  a  factor 
of  25  when  the  distance  increased  from  1  to  7  cm;  for  experiments  #8-1 
and  11-1  at  a  pressure  of  500  microns  there  is  a  decrease  of  "400  in 
going  1.7"  (3.8  cm)  from  the  plasma.  This  suggests  the  importance  of 
wall  collision  neutralization.  When  a  short  set  of  copper  probes  were 
used  (the  tips  did  not  penetrate  into  the  discharge  tube)  a 
characteristic  curve  was  not  obtained;  presumably  the  ion  density  was 
below  the  critical  limit  to  form  a  sheath  around  the  probes.  The 
implication  is  that  the  restriction  at  the  end  of  the  discharge  tube 
will  drastically  reduce  the  number  of  ions  that  can  flow  into  the 
reaction  vessel. 

The  ion  densities  increase  with  increasing  microwave  power 
(experiments  #1-1  and  2-1  or  4-1,  5-1,  and  6-1,  or  7-1  and  8-1  or  10-1 
thru  14-1);  however,  the  electron  temperature  of  the  hot  electrons 
increase  with  power  up  to  "80%  and  then  decreases  as  the  power  is 
increased  to  100%  (experiments  10-1  thru  13-1  and  14-1).  This  decrease 
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can  be  explained  by  the  loss  of  available  power  to  the  discharge  tube 
walls.  Recently  it  has  been  reported34  that  when  the  plasma  plume  is  in 
direct  contact  with  the  wall,  the  wall  erodes.  The  energy  necessary  for 
this  errosion  (decomposition)  is  provided  by  the  microwave  generator, 
thus  decreasing  the  energy  available  to  produce  hot  electrons. 

Lower  electron  temperatures  and  densities  are  the  apparent  effect 
of  using  larger  diameter  probes  as  was  observed  by  replacing  the  smaller 
copper  probes  with  the  larger  tungsten  set.  This  has  also  been  reported 
by  Cohen35  and  Little  and  Waymouth3  .  Cohen35  has  shown  that  the  ratio 
of  the  radius  of  the  wall  boundary  (Ry)  to  the  probe  radius  (Rp) 
determines  the  electron  temperature  and  density.  For  Rw/Rp=100  the 
calculated  temperatures  are  20%  less  than  for  an  infinite  ratio  while 
the  densities  vary  by  only  1-2%;  Rw/Rp=  50  will  produce  a  50%  decrease 
in  the  calculated  temperature.  In  our ^experiments  the  tungsten  probes 
have  Rw/Rp  <  50  so  it  is  expected  that  the  calculated  temperatures  could 
be  less  than  50%  of  the  true  temperature.  The  shorter  penetration  of 
the  tungsten  probe  into  the  plasma  could  also  account  for  the  smaller 
electron  density. 

Reaction  Vessel  Results; 

Our  earlier  work3®  indicated  ion  concentrations  <0.03%  of  the 
atom  contentration;  this  is  based  on  the  assumption  that  the 
density  of  negative  ions  is  equal  to  the  density  of  all  available 
electrons,  i.e.  an  upper  limit  for  the  number  of  negative  ions  is 
determined.  By  assuming  50%  dissociation  of  the  H2  then  the  H  atom  flow 
through  the  exit  orifice  of  the  discharge  tube  is  equal  to  the  measured 
H2  flow  into  the  discharge  tube.  Having  shown  the  agreement  of  the 
double  floating  probe  technique  in  the  cavity  region  the  present  work 
also  indicates  that  ±he  negative  ion  concentration  is  less  than  1.8  * 
10-4  %  and  9.8  *  10-5  %  for  96  and  60%  power,  respectively. 

Additionally,  the  electron  temperature  measured  in  the  reaction  vessel 
are  higher  than  in  the  discharge  tube;  this  may  be  due  to  the  radius 
ratio  of  “500  as  compared  to  50. 

RESIDENCE  TIMES;  Pressure  time  profiles 


^erimental  Setup; 


The  diatomic  reactant  is  pulsed  using  a  Bosch  fuel  injector 
(previously34  a  homemade  solenoid  valve  was  used);  the  on/off  duty 
cycle  and  cycletimes  can  be  varied.  A  nude  ionization  gauge  (NRC  #538) 
was  used  for  measuring  the  pressure-time  profiles.  The  gauge  was 
mounted  on  a  glass  rod  so  that  the  position  inside  the  reaction  vessel 
could  be  changed.  A  glass  encapsulated  ionization  gauge  was  mounted  on 
top  of  the  reaction  gauge;  its  conductance  was  limited  by  a  1/2" 
diameter  entrance  tube  2"  long. 


A  standard  ionization  gauge  power  supply  was  used  for  biasing  the 
grid  and  collector  and  to  provide  current  to  the  heater  for  degassing. 
The  collector  was  connected  to  an  operational  amplifier  configured  as  a 
current  to  voltage  converter.  This  output  was  then  input  to  a  variable 
gain  inverting  operational  amplifier  and  subsequently  fed  to  an 
oscilloscope  or  PARC/EG&G  model  4203  signal  averager.  The  signal 
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averager  was  used  to  reduce  any  60  Hz  noise.  The  signal  was  averaged 
for  25-60  pulses:  2048  data  points  were  used  to  represent  the  time 
profile.  The  time  response  of  the  gauge  and  amplifiers  was  <0.1  msec. 
The  data  from  the  signal  averager  was  then  transferred  to  a 
microcomputer  for  data  formatting  and  then  to  a  Prime  750  mainframe 
computer  for  data  processing  and  plotting. 

Results  and  Discussion: 

The  pressure-time  profiles  were  analyzed  with  a  non-linear  least 
squares  program;  the  starting  values  for  this  program  were  estimated  by 
a  zero  order  linear  least  squares  analysis.  The  general  profile 
indicated  three  regions:  a  rise,  a  fast  exponential  decay  and  a  slower 
exponential  decay.  The  pumpout  rate  constants  (i.e.  the  slower  decay 
rates)  were  determined  using  data  from  point  number  1000  to  2000.  In 
general,  the  peak  heights  increase  with  flow;  the  response  for  argon  is 
greater  than  that  of  hydrogen  due  to  the  larger  ionization  efficiency 
for  argon. 

Two  types  of  experiments  were  done:  with  and  without  the  inner  can 
(normally  in  a  CM  experiment  the  inner  can  is  cooled  with  liquid 
nitrogen  to  facilitate  complete  vibrational  deactivation  of  the  excited 
reaction  product  which  collides  with  the  wall).  Experiments  with  argon 
and  hydrogen  were  performed  for  various  duty  cycles  and  cycle  times. 

The  nude  gauge  was  placed  at  either  the  center  of  the  vessel  or  at  the 
end  while  the  encapsulated  gauge  was  permanently  located  at  the  top 
(near  center)  of  the  reaction  vessel.  The  "on"  time  (time  the  valve  was 
opened)  varied  from  10  to  20  milliseconds  and  the  cycle  time  (the  time 
between  pulses)  was  either  100  or  200  milliseconds. 

Figure  II-3  illustrates  pressure  time  profiles  for  hydrogen  with 
the  inner  can  and  "on"  times  of  10,  15  and  20  msecs  (the  total  flow  rate 
remained  constant).  The  starting  trigger  times  are  the  same;  the  delay 
between  the  energiziation  of  the  injector  and  the  detection  of  the 
pressure  pulse  is  due  to  the  mechanical  delay  of  the  injector.  The  peak 
times  are  equal  to  the  time  of  the  energization.  When  the  cycle  time  is 
increased  to  200  msec  (see  Figure  II-4)  the  baseline  approaches  zero, 
indicating  that  the  vessel  has  been  completely  pumped  out.  With  the 
inner  can  in  places  double  exponential  for  the  decay  portion  of  the 
profile  is  observed.  This  is  somewhat  masked  when  the  can  is  removed  as 
shown  in  Figure  II-4  for  both  argon  and  hydrogen. 

When  the  cycle  time  is  decreased  to  100  msec  the  baseline  increases 
and  a  decay  between  the  trigger  and  before  the  next  pressure  increase  is 
observed.  When  the  nude  gauge  is  placed  at  the  end  of  the  vessel  the 
pressure  profile  exhibits  a  rounded  top;  due  to  the  longer  transit  time 
the  peak  is  broadened  at  the  top.  When  the  encapsulated  gauge  (which 
has  an  inherently  small  pumping  conductance)  is  used.  Presumably  the 
time  constant  of  this  gauge  has  been  substantially  increased  and  the 
gauge  container  is  acting  as  a  ballast.  The  experimental  conditions  and 
calculated  results  for  all  of  the  pressure  profile  experiments  are 
tabulated  in  Table  II-2. 

In  general  the  pumpo.pt  out  rate  for  H2  (20  sec-1)  is  larger  than 
that  for  argon  (  “16  sec-1).  The  differences  between  the  pump  out  rate 
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Figure  1 1 -3 .  Pressure-time  profiles  with  nude  gauge  located  below 
pulser  for  hydrogen  with  the  inner  can  and  a  cycle  time  of  100 
msec.  The  "on"  times  are:  10  msec  (curve  1,  experiment  *1-P),  15 
msec  (curve  2,  experiment  #2P),  20  msec  (curve  3,  experiment  *3- 
P). 
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Figure  1 1 -4 .  Pressure  time  profiles  with  nude  gauge  located  below 

pulser  with  a  cycle  time  of  200  msec  and  "on"  times  of  10  msec. 

curve  1:  hydrogen  with  inner  can,  experiment  *4-P 

curve  2:  hydrogen  without  inner  can,  experiment  «5-P 

curve  3:  argon  with  inner  can,  experiment  #6-P 

curve  4:  argon  without  inner  can,  experiment  S7-P 


constants  with  and  without  the  inner  can  are  small  and  within 
experimental  error/  thus  they  are  not  considered  distinguishable. 
However  the  double  exponential  decay  is  more  apparent  with  the  inner 
can.  Presumably  this  double  exponential  decay  can  be  attributed  to  the 
possibility  that  molecules  which  leave  the  inner  can  return  before 
eventually  exiting  through  the  diffusion  pump. 


CONCLUSIONS 

These  diagnostic  experiments  indicate  that  the  ions  produced  from  a 
microwave  discharge  and  entering  a  reaction  vessel  through  a  restriction 
are  not  at  a  sufficiently  large  concentration  to  allow  NIMR's  to 
interfere  with  atomic  reactions.  The  relative  flows  of  ions  was 
determined  to  be  less  than  1.8  10“°  times  the  flow  of  atoms,  so  that 
even  with  a  factor  of  1000  increase  in  reaction  cross  section  NIMR's 
will  not  be  competitive  with  the  atom-diatomic  reactions. 

The  pressureprof iles  indicate  that  the  pumping  speed  for  this 
system  is  “20  sec-  ,  giving  a  residence  time  of  50  msec.  For  the  CM 
experiments  radiative  relaxation  for  HF  is  important  on  this  time  scale; 
for  steady  state  experiments  a  background  'soup'  comprised  of  products 
and  reactants  will  build  up  so  that  secondary  reactions  and  collisional 
relaxation  will  be  important.  although  these  diagnostics  were 

performed  for  a  specific  CM  experimental  apparatus  the  results  are 
applicable  to  other  experiments  in  which  atoms  are  formed  in  a  microwave 
dischare  and/or  reactants  are  pulsed  in  a  low  pressure  environment. 
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III.  CM  MODELING 


INTRODUCTION 

Chemiluminescent  Mapping  (CM)  has  been  developed  as  an  alternative 
method  in  obtaining  formation  and  relaxation  information  for  exoergic 
bimolecular  exchange  reactions.  Both  time  and  spectral  resolution  are 
obtained  simultaneously.  Conventionally  steady  state  experiments  have 
been  performed  (Arrested  and  Measured  Relaxation  and  Fast  Flow 
experiments)  in  which  constant  flows  of  the  atomic  and  molecular 
reactants  are  maintained  and  the  chemiluminescence  is  recorded.  A 
typical  reaction  is 

F  +  H2  - >  HFV/J  +  H  kVrJ  (v-1-3) 

Two  disadvantages  of  the  steady  state  experiments  are:  i)  that  for  slow 
reactions  the  gas  flows  must  be  increased  to  obtain  a  useable  emission 
intensity  which  in  turn  requires  larger  pumping  speeds  to  maintain  the 
required  low  pressure,  and  ii)  that  the  observed  luminescence  signal  is 
a  composite  of  blackbody  background  and  the  desired  chemiluminescence 
emission.  For  reactions  that  have  large  exoergicities  the  scanned 
spectral  region  is  sufficiently  large  so  that  the  change  in  blackbody 
background  can  be  greater  than  or  comparable  to  the  emission  signal  for 
the  reaction.  However,  compensation  for  the  background  can  be  made  by 
internally  modulating  the  signal,  i.e.  by  pulsing  either  one  of  the 
reactants.  Variations  of  the  pulse  rate  and  duty  cycle  allow  collection 
of  different  types  of  information.  A  slow  pulse  rate  and  a  long 
reaction  "on"  time  is  equivalent  to  two  experiments.  With  the  reactant 
"on"  the  signal  plus  background  is  detected,  while  with  the  reactant 
"off"  only  the  background  is  detected.  The  numerical  difference  between 
these  signals  is  that  which  is  due  solely  to  the  emission  produced  by 
the  reaction.  Extensive  signal  averaging  beyond  that  of  phase  sensitive 
detection  is  easily  implemented  for  CM  experiments. 

A  single  CM  pulse  cycle  consists  of  two  stages:  input  (reactant 
injection)  and  output  (relaxation  of  reactants  and  or  products).  Two 
limiting  regions  can  be  considered  for  the  input  cycle:  i)  instantaneous 
input  in  which  the  formation  time  is  substantially  less  than  relaxation 
processes  (i.e.  reaction)  and  ii)  slow  input  in  which  the  formation  time 
is  slow  relative  to  reaction  times.  Further  simplifications  can  be  made 
by  requiring  that  one  of  the  reactants  is  in  excess  (this  reduces  second 
order  kinetics  to  first  order  kinetics).  Thus  during  a  CM  cycle  one  of 
the  reactants  is  depleted  via  reaction  while  the  excited  product  is 
depleted  via  radiative  or  collisional  relaxation,  physical  removal  by 
pump  out  or  deactivated  by  wall  collisions  to  its  ground  state.  The 
relative  importance  of  each  of  the  above  processes  is  related  to  the 
initial  conditions  and  reactant  preparation  time. 

A  CM  experiment  can  be  performed  by  pulsing  either  the  atomic  or 
diatomic  (molecular)  species;  both  techniques  have  specific  advantages. 
Our  earlier  CM  experiments  were  based  on  pulsing  the  diatomic  molecule. 
In  this  case  there  is  an  initial  increase  in  pressure  followed  by  a 
decrease  for  each  cyle.  In  this  report  we  denote  this  CM  Pressure  Pul se 
experiment  as  CM-PP.  Alternatively  in  a  CM-LP  (iaser  Pulse)  experiment 
the  atomic  reactant  is  pulsed;  using  laser  excitation  the  atom  is  formed 
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nearly  instantaneously  under  approximate  isobaric  conditions. 

Four  different  time  scales  must  be  considered  in  a  CM  experiment: 
reactant  formation,  reaction,  relaxation,  and  sampling  (observation). 

It  is  more  difficult  to  obtain  information  about  a  specific  process  if 
the  time  scales  are  compar  abl  e  than  it  is  if  the  time  scales  are  well 
seperated.  It  is  also  important  that  the  sampling  interval  time  is 
short  enough  to  resolve  the  fastest  process  otherwise  significant  time 
averaging  will  result.  If  sufficient  relaxation  occurs  during  a  single 
CM  cycle,  the  relative  composition  of  the  reaction  mixture  can  change 
dramatically.  This  can  then  cause  a  change  in  the  relative  importance 
of  the  different  relaxation  processes.  Computer  modeling  aids  in  the 
understanding  of  which  processes  are  dominant  in  the  different  time 
regions  of  a  CM  cycle.  This  knowledge  can  be  used  to  optimize  the 
experimental  conditions  necessary  to  obtain  specific  information 
(nascent  energy  distribution,  total  reaction  rate  constants, 
deactivation  rate  constants). 

It  is  also  important  that  the  ensemble  of  species  which  is  being 
observed  is  homogeneous.  This  corresponds  to  a  uniform  concentration 
over  the  observation  volume.  The  homogeneity  can  be  obtained  if 
reaction  takes  place  in  a  small  volume.  The  observed  emission  signal 
which  is  directly  proportional  to  the  observation  volume  is  degraded 
resulting  in  a  lower  signal  to  noise  ratio.  Thus  a  trade  off  between 
small  observation  volume  and  emission  intensity  must  be  made.  The 
calculations  presented  here  assume  a  homogeneous  mixture  so  that  spatial 
averaging  is  not  necessary. 

The  first  generation  CM  exper iments14/1®,26a  were  based  on  pulsing 
the  molecular  reagent  (CM-PP).  Initially,  we  modeled  these  experiments 
and  the  simulation  results  indicated  that  experimental  modifications 
were  needed.  Simulations  that  incorporated  these  modifications  into  the 
CM-PP  model  were  then  performed.  Additional  model  computer  simulations 
were  then  performed  for  both  CM-PP  and  CM-LP  experiments;  these  results 
were  then  used  to  fine  tune  the  experiments  so  that  nascent  energy 
distributions,  total  rate  constants  and  relaxation  rate  constants  could 
be  extracted  from  the  experimental  data.  The  results  of  these  CM-PP  and 
CM-LP  model  calculations  are  discussed  separately  in  this  section.  The 
'new'  experiments  are  reported  in  section  V 

CM-PP  MODEL  CALCULATIONS 


The  CM-PP  experiments  consist  of  F  atoms  (produced  in  a  microwave 
discharge  of  CF^)  with  a  flow  rate  F^a  into  a  low  pressure  reaction 
vessel  with  liquid  nitrogen  cooled  walls.  The  reaction  vessel  is 
connected  to  a  diffusion  pump  with  a  first  order  pumping  rate  constant, 
kp.  Hydrogen  is  pulsed  into  the  reaction  vessel  at  a  rate,  H2^n,  with 
an  "on"  time  of  tpp  and  a  total  cycle  time  of  tcyc.  The  F  atoms  can 
react  with  H2  to  form  HF(v)  (for  simplicity  the  rotational  quantum 
number,  J,  has  been  dropped),  diffuse  to  the  walls  of  the  vessel  to  form 
F or  be  physically  removed  from  the  reaction  vessel  (i.e.  pumped  out). 
Likewise  H2  can  react  with  F  atoms  to  form  HF(v),  deactivate  HF(v),  or 
be  pumped  out  of  the  vessel.  The  excited  product,  HF(v)  can  be  depleted 
by  radiative  or  collisional  relaxation,  deactivated  by  wall  collisions 
to  its  ground  state,  or  be  pumped  out  of  the  vessel.  The  following 
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reaction  scheme  is  used  to  model  the  CM-PP  experiments: 


ki 

F  +  H2  - >  HF  ( i)  +  H 

kp 

F  - e - > 

k^ 


HF  ( i) 

*w 

- > 

HF  (0) 

HF  ( i) 

+  h2 

kj,i<H2> 

HF  ( j ) 

+  h2 

HF  (i) 

+  HF(0) 

HF  ( j ) 

+  HF  (k ) 

HF  ( i) 

HF  ( j ) 

This  reaction  scheme  can  be  described  by  the  following  set  of 
differential  equations  in  which  it  is  assumed  that  i  =  j  +/-  1  (only 
nearest  neighbor  transitions  are  allowed). 

d[F]/dt  = 

Fin  -  ktotfH2nF]  -  kp[F]  -  kw[F] 
d[H2]/dt  =  H2in  -  ktot [H2]  [F]  -  kp[H2] 
d  [HF  (i)  ]  /  dt  =  ki[H2][F]  +  Airj[KF(j)] 

+  kif  j  (H2)  [H2]  (HF  ( j  )  ]  -  Aj/jL[HF(i)] 

-  kjfi(H2)  tH2l(HF(i)]  -  kplHF(i)] 

-kw  [HF  (  i)  ]  +  kifj(HF0)  [HF  ( j  )  ]  [HF(0)  ] 
-kjri(HF0)  [HF ( i) )  IHF{0) ] 

=  microscopic  reaction  rate  constant  for  state  i 
ktot  “  ki 

Aj  ^  =  Einstein  coefficient  for  emission  from  state  j 
'  to  state  i 

kff  j(H2)  =  HF-H2  bimolecular  collisional  energy 
? "  transfer  rate  constant  from  state  j  to  i 

k  (HFq)  =  HF-HF(O)  bimolecular  collisional  energy 
f  transfer  rate  constant  from  state  j  to  i 

kp  =  pump  out  rate  constant 
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kw  =  wall  collision  rate  constant 


These  differential  equations  are  numerically  integrated  to  obtain 
[H2],  [F],  and  [HF (v)  ]  as  a  function  of  time  for  various  experimental 
conditions  and  rate  constants.  An  IMSL  (International  Mathematical  and 
Statistical  Libraries)  subroutine  (Runge-Kutta-Verner  5th  and  6th  order 
Method)  is  used  for  the  numerical  integration;  the  convergence  tolerance 
was  set  for  0.001  while  the  time  interval  ranged  from  5  x  10“®  to  1  x 
10-4  sec  depending  on  the  cycle  time. 

Two  experimental  variables  are  F^p  and  H2j_.  The  absolute  [F]  is 
unknown,  but  will  be  less  than  four  times  the  CF^  flow  rate  divided  by 
kp  +  k„.  Hydrogen  pressure  profiles  have  been  measured2®13  and  are 
character ized  by  a  fast  rise,  initial  fast  exponential  decay,  and  a 
second  slower  exponential  decay.  A  pump  out  rate  constant,  kp,  was 
determined  by  a  linear  least  squares  fit  of  the  slower  exponential 
decay2  .  To  be  noted  is  that  there  are  equal  pressures  of  H ?  during 
the  input  ("up")  and  output  ("down")  stage.  Differences  in  the  HF(v) 
emission  intensity  at  these  equivalent  pressure  points  can  be  attributed 
to  changes  in  the  relative  composition  of  the  environment. 

The  H2  pressure-time  profile  was  simulated  using  a  double 
exponential  function  with  tQn  =  10  msec  and  tcyc  =  100  msec.  A  linear 
least  squares  fit  of  the  slow  exponential  decay  was  used  to  obtain  a  k 
and  the  "best"  simulated  H2  pressure-time  profile  was  chosen  to  be  they 
one  in  which  the  kp  obtained  was  equivalent  to  the  experimental  kp. 
Either  the  experimental  or  parameterized  H2  pressure-time  profiles  were 
used  in  solving  the  above  differential  equations  with  no  meaningful 
difference  in  the  results. 

Typical  HF  emission  intensity-time  profiles  obtained  from  diagnostic 
CM  experiments39  are  shown  in  Figure  III-l.  These  profiles  are  for 
moderate  (10  umoles/sec)  and  high  (40  umoles/sec)  flows  of  H2  with  a  tpn 
=  10  msec  and  tcyc  =  100  msec.  At  moderate  flows,  the  HF  follows  the 
hydrogen  profiled  peaking  when  the  H2  is  turned  off.  At  high  flows  of 
H2,  the  HF  intensity-time  profiles  exhibit  a  double  maxima.  The  first 
maximum  and  decay  occurs  while  the  H2  is  "on".  After  the  H2  is  turned 
off,  there  is  another  increase  and  final  decay. 

In  Figure  III-2,  the  data  from  Figure  III-l  is  plotted  as  a  function 
of  the  experimentally  measured  Ho  pressure.  This  type  of  graph 
indicates  how  HF(v)  varies  with  the  H2  environment.  Because  there  are 
equivalent  H2  pressures  during  the  input  and  pump  out  stages,  it  is 
expected  that  the  "up"  and  "down"  branches  of  this  plot  will  be 
super imposabl e  for  equivalent  amounts  of  available  H2.  However,  the 
experimentally  observed  time  dependence  of  H2  is  experimentally  measured 
without  the  presence  of  F  so  that  the  experimental  H2  profile  is  for  the 
total  H2  input  to  the  system  and  not  the  amount  of  available  H2.  For 
small  percent  reaction,  the  difference  between  these  quantities  is 
negligible.  However,  this  difference  will  account  for  a  slight 
displacement  of  the  "up"  and  "down"  branches  (in  general,  the  "down" 
branch  is  displaced  to  the  right  of  the  "up"  branch).  In  the  plots  of 
Figure  III-2  both  a  displacement  between  the  "up"  and  "down"  branches 
and  a  change  of  curvature  are  to  be  noted. 


Figure  III-l.  Experimental  HF(v)  Emission  Intensity  Versus 
Time,  t  =10  msec,  t  =  100  msec 


Parameters  in  the  simulation  were  varied  in  order  to  qualitatively 


reproduce  the  features  of  these  experimentally  observed  profiles.  The 
parameters  used  were  Ag  ■,  =  180  sec”1,  Ai  2  =  320  sec”1,  A2  3  =  400 
sec”1,  ktot=  1  x  10°  torr  1  sec,  k3/k2  =t.5,  ki/k2  =  0.3,  'k  •  ^(H2)  = 

2  x  105  torr”1  sec”1,  k_  =  40  sec”1,  kw  d  ■,  =  500  sec  -1,  k-,  2(HF0^  =  5.0 
x  105  torr  1  sec”1,  k2  3(HFg)  =  9.0  x  llT'torr”1  sec”1.  Tne  Einstein 
coefficients  are  from  reference  4,  k^0^  from  reference  40,  k3/k2  and 
ki/k2  from  reference  7,  and  k^j(HFg)  from  reference  41. 

Simulation  Results  And  Discussion 

Simulations  of  the  first  generation  CM-PP  experiments  (tQn  =  10 
msec)  and  the  modified  CM-PP  experiments  (tQn  =  0.5msec)  are  discussed 
in  this  section.  The  discussion  is  divided  into  four  parts.  Each  part 
is  identified  by  a  number:  "1"  (low  reagent  flow)  or  "2"  (high  reagent 
flow)  or  a  letter:  "a"  (tQn  =  10  msec)  or  "b"  (tQn  =  0.5  msec) 

la.  Low  Flow  (F/H2  >  30),  tQn  =  10  msec 

In  Figure  III-3,  the  time  dependence  of  the  partial  pressures  of  F 
and  H2  is  displayed.  The  F  atoms  are  17%  depleted  and  H2  is  the 
limiting  reagent.  The  H2  partial  pressure  initially  increases  while 
that  for  F  decreases.  When  the  H2  is  turned  off,  the  H2  decays  while 
the  F  increases  back  to  its  non-reactive  steady  state  value.  The  H2  is 
depleted  by  reaction  and  pumping  and  this  depletion  increases  with  time. 

The  HF (v)  pressure-time  profiles  are  shown  in  Figure  III-4  and  mimic 
the  H2  pressure  time  profiles.  Also,  for  this  cycle  time,  significant 
background  is  present  (the  curves  do  not  relax  to  zero  pressure) 
indicating  that  the  reaction  vessel  is  not  completely  pumped  out  before 
the  next  H2  pulse. 

A  comparision  between  equivalent  steady  state  experiments  and  a 
single  CM-PP  experiment  can  be  made  for  the  time  dependence  of  [F]  or 
the  dependence  of  [F]  with  H2  pressure.  The  deviation  between  the  "up" 
and  "down"  branches  of  the  steady  state  calculations  is  due  to  the  fact 
that  the  total  H2  pressure  is  used  and  not  the  instantaneous  H2  pressure 
which  takes  into  account  the  H2  that  has  reacted.  If  there  is  a  small 
percent  reaction,  the  deviation  between  the  two  branches  is  reduced.  The 
deviation  for  the  real  time  dependence  of  F  with  H2  pressure  is  larger 
than  that  observed  for  the  steady  state  curves.  Tne  finite  time  for 
reaction  and  relaxation  accounts  for  this  difference.  Similar 
dependences  for  the  HF(v)  partial  pressures  are  plotted  in  Figure  III-5 
for  HF  ( v  =  3) . 

Molecular  and  atomic  fluxes  are  the  individual  rate  terms  in  the 
differential  equations  for  H2,  F  and  HF(i)  and  their  magnitudes  indicate 
which  processes  are  dominant  at  given  times  during  a  CM  cycle.  Tables 
III-l  thru  III-3  summarize  these  fluxes  for  HF(v)  at  low  and  high  flew 
conditions.  For  low  flows,  the  input  (production)  terms  dominate  so 
that  the  net  flux  is  positive  until  the  H2  is  turned  off,  then  the 
output  (depletion)  terms  (net  flux  negative)  cause  the  HF  to  decay.  The 
dominant  processes  are  reaction,  collisions  with  the  vessel  walls,  and 
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HFty=3)  PRESSURE  (torr) 


Figure 
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1 1 1 - 5 .  HF(v=3)  Pressure  Versus  H«  Pressure, 
CM-PP  Simulations:  Low  Flow,  t  =10 
t _  =  100  msec  on 
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radiative  relaxation.  Collisions  with  H2  and  HF(v=0)  are  negligible  for 
low  flow  conditions. 

Fluorine  atoms  were  found  to  be  depleted  mainly  by  recombination  on 
the  walls  of  the  vessel  followed  by  reaction  and  pump  out  from  the 
vessel.  Wall  collisions  account  for  “54%  of  HF(v=3)  depletion  followed 
by  radiative  relaxation  (“42%).  Because  the  Einstein  coefficients 
become  smaller  with  decreasing  vibrational  level,  the  fraction  of 
depletion  due  to  radiative  relaxation  decreases  as  the  vibrational 
energy  is  decreased.  Wall  collisions  account  for  “58%  of  the  depletion 
of  HF(v=2)  while  “37%  is  due  to  radiative  emission.  For  HF(v=l),  “69% 
of  the  depletion  is  by  collisions  with  the  walls  and  “25%  is  due  to 
radiative  relaxation. 

Reaction  is  the  only  input  to  HF(v=3)  and  is  the  dominant  input 
(80%)  to  HF(v=2).  Radiative  relaxation  from  HF(v=3)  accounts  for  only 
20%  of  the  input  to  HF(v=2).  Initially  radiative  relaxation  from 
HF(v=2)  into  HF(v=l)  is  dominant,  between  2  and  5  msec  reaction  becomes 
more  important,  and  finally  radiative  relaxation  is  again  dominant.  If 
the  cycle  time  is  increased,  reaction  is  the  dominant  input  to  HF(v=l) 
until  5  msec  when  radiative  relaxation  from  HF(v=2)  takes  over.  These 
results,  again,  illustrate  that  the  reaction  vessel  is  not  completely 
evacuated  between  pulses  for  tcyc  =  100  msec.  For  radiative  relaxation 
to  be  the  initial  dominant  input,  residual  HF(v=2)  must  be  present. 
However,  even  when  the  cycle  time  is  increased,  radiative  relaxation  is 
a  dominant  input  to  HF(v=l)  after  5  msec  in  contrast  to  HF(v=2)  where 
reaction  is  dominant  for  all  times.  This  is  due  to  the  fact  that  the 
microscopic  reaction  rate  constant  for  HF(v=2)  is  larger  than  for 
HF  ( v  =  l)  or  HF(v  =  3)  (k^kj.  =  3.3,  k2/k3  =  2). 

Figure  III-6  shows  plots  of  HF (v)/HF ( v=2)  as  a  function  of  time  for 
both  tcyc  =  100  msec  and  200  msec.  For  a  cycle  time  of  200  msec,  the 
HF (v)  and  H2  decay  to  zero  and  the  ratio  plots  are  linear  out  to  2  msec. 
An  extrapolation  of  the  linear  portion  of  these  curves  to  zero  time 
gives  the  nascent  distributions  (k 3/k 2  =  0.5,  k^/ko  =  0.3).  The  slopes 
of  the  curves  depend  on  the  relaxation  processes;  k 3/k 2  has  a  negative 
slope  and  k^/k2  has  a  positive  slope.  Unfortunately,  only  data  for  20% 
of  tQQ  can  be  used  for  the  extrapolation  so  that  the  majority  of  the 
data  does  not  provide  direct  information  on  the  nascent  distributions. 
These  initial  points  will  have  the  lowest  intensity  and  the  smallest 
signal/noise  ratio  so  that  the  extrapolated  distribution  would  have  very 
large  errors  associated  with  it.  The  ratio  plots  obtained  for  a  cycle 
time  of  100  msec  indicate  that  if  the  vessel  is  not  completely  evacuated 
between  pulses,  nascent  distributions  are  not  easily  obtained. 


2a.  High  Flows  (F/H2  <  0.75),  tQn  =  10  msec 

The  time  dependence  of  the  partial  pressures  of  F  and  H2  illustrate 
that  F  atoms  are  the  limiting  reagent  (“80%  depletion).  Reaction  is  now 
the  dominant  depletion  process  of  F  at  the  expense  of  wall 
recombination.  The  difference  between  the  H2  profiles  with  and  without 
reaction  is  small  due  to  the  small  percent  reaction.  Figure  IV-7  shows 
the  time  dependence  of  the  HF(v).  Double  maxima  are  observed;  the  first 
appears  before  the  H2  is  turned  off  and  the  second  after. 
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The  [F]  variation  with  H2  pressure  is  similar  to  that  observed  for 
low  flows.  The  deviation  between  the  "up"  and  "down"  curves  for  the 
steady  state  calculations  of  [F]  is  smaller  because  of  less  percent 
reaction.  Figure  IV-8  shows  the  HF(v=3)-H2  profile  which  is  also 
similar  to  the  behavior  of  HF(v=l)  and  HF(v=2).  Steady  state 
calculations  exhibit  a  maximum  before  and  after  the  H2  is  turned  off. 
This  behavior  can  be  understood  by  considering  the  steady  state  equation 
for  HF(v=3): 


HFSS  (v=3)  =  { (k3  tH2]  )/(k2jr3  IH2]  +  kr)}{Fin/(ktot^H25  +  kpH  HI-1) 

where  kr  is  the  sum  of  all  depletion  processes  that  are  independent  of 
H2*  Because  [H2]  varies  with  time,  the  importance  of  H2  collisional 
deactivation  (k2  3  )  and  F  atom  depletion  (kx.pt[H2])  will  also  vary 

with  time.  For  'small  [f^l  {early  times),  k2  3IH3]  <<  kr  (absence  of  Ho 
collisional  deactivation)  and  <<  k-'Tminimal  F  atom  depletion), 

so  that  the  increase  of  HFss(v=3)  is  proportional  to  [H2 ] -  For  large 
[H2],  k?  3  [H2 ]  >>  k  and  ktot[H23  >>  kp  so  that  HFss(v=3)  is 
proporti6nal  to  1/[h2]»  Since  IH2]  is^large  and  increasing  for  times  < 
tDn,  the  HFss(v=3)  decreases.  The  maximum  in  Figure  III-8  occurs  when 
the  denominator  is  proportional  to  tH23;  two  possibilities  are  either 
k2  3  [ H2 ]  >>  k  or  ktotJH2^  >:>  kp  so  that  HFss(v=3)  is  independent  of 
[H21*  The  behavior  or  tne  "dowh"  branch  in  Figure  III-8  is  similar  to 
the  "up"  branch  since  there  are  equivalent  H2  pressures  during  the  input 
and  pump  out  stages.  The  "real"  HF(v)  builds  up  to  values  greater  than 
those  for  a  steady  state  due  to  insufficient  time  for  the  F  atom 
concentration  to  be  depleted  and  collisional  relaxation  of  HF(v). 


Collisional  deactivation  of  HF(v)  by  H2  is  of  comparable  magnitude 
to  radiative  relaxation  and  collisions  with  the  walls  (see  Tables  III-4 
thru  III-6).  Initially  collisions  of  HF(v=3)  with  the  wall  dominate 
over  depletion  by  radiative  or  collisional  relaxation  as  observed  for 
low  flows.  While  the  H2  is  "on"  collisional  deactivation  of  HF(v)  by 
H2  becomes  dominant  and  accounts  for  up  to  50%  of  the  relaxation. 
Similar  results  for  the  depletion  of  HF(v=l)  and  HF(v=2)  are  observed. 


Input  to  HF(v=2)  is  dominated  by  reaction  (~80%)  at  all  times  as  it 
was  for  the  low  flow  conditions;  input  due  to  H2  collisional  relaxation 
from  HF(v=3)  is  only  ~20%.  As  in  the  low  flow  case,  the  effect  of 
residual  HF(v)  is  illustrated  by  the  fact  that  radiative  relaxation  from 
HF(v=2)  is  the  initial  input  to  HF(v=l).  If  the  cycle  time  is 
increased,  reaction  is  the  initial  dominant  input.  The  main  input  to 
HF(v=l)  after  approximately  3  msec  is  collisional  relaxation  from 
HF(v=2). 


The  HF(v)/HF(v=2)  ratio  as  a  function  of  time  is  shown  in  Figure 
III-9.  As  was  observed  for  low  flow  conditions,  the  ratio  plots  for 
tgy„  =  100  msec  cannot  be  used  to  obtain  nascent  distributions  because 
ofythe  presence  of  residual  HF(v).  If  the  cycle  time  is  increased  to 
200  msec,  the  plots  show  linear  behavior  during  the  time  when  relaxation 
processes  are  of  less  importance  than  reaction  (16%  of  ton).  Again  the 
major  portion  of  the  data  does  not  provide  information  that  can  be 
readily  interpreted. 
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The  results  of  the  above  two  sets  of  simulations  indicate  that  in 
order  to  obtain  nascent  distibutions,  the  cycle  time  must  be  long  enough 
to  ensure  that  there  is  no  residual  HF(v)  present  and  that  relaxation 
processes  must  be  minimized  during  the  observation  time.  The  latter 
suggests  that  a  shorter  pressure  pulse  and  usable  data  for  a  given 
observation  time  is  maximized.  Therefore,  the  "old"  CM-PP  experimental 
system  needed  to  be  modified;  a  second  set  of  simulations  were  done  to 
show  that  the  "new"  experiments  would  be  an  improvement.  The  "new" 
changes  in  the  model  are  that  the  H2  is  pulsed  into  the  reaction  vessel 
with  a  tQn  =  0.5  msec  and  a  tc.,c  =  TO  msec.  In  addition,  the 
observation  volume  was  reduced  so  that  HF(v)  is  homogeneously  produced 


in  this  volume;  reagents  flow  out  of  this  volume  with  a  rate  constant, 
kp.  Because  of  the  smaller  observation  time  and  the  reduced  observation 
vfelume,  HF(v)  does  not  have  time  to  collide  with  the  vessel  walls  and 
return  to  the  observation  region;  the  kv  term  was  omitted.  The  value  of 
k  used  for  HF(v)  flow  from  the  observation  volume  was  1.3  x  103  sec-1 
while  all  other  parameters  remained  the  same. 


lb.  Low  Flows  (F/H2  >  30),  tQn  =  0.5  msec 

The  Ho  pressure  increases  with  time  and  at  0.5  msec  starts  to 
decrease.  There  is  less  percent  reaction  than  was  observed  for  the  same 
flow  conditions  and  tQn  =  10  msec.  The  F  atoms  decay  (  25%  depletion) 
but  do  not  reach  a  minimum  until  after  the  hydrogen  is  turned  off. 
Depletion  of  F  atoms  is  due  mainly  to  wall  recombination  followed  by 
reaction  and  pump  out  of  the  reaction  vessel.  The  HF(v)  partial 
pressure  -  time  profiles  show  a  single  maximum  that  does  not  occur  until 
after  the  H2  is  turned  off. 

The  HF(v)-H2  and  F-H2  pressure  profiles  show  a  greater  deviation 
from  the  steady  state  calculations  than  was  observed  for  the  ton  =  10 
msec  case.  This  is  to  be  expected  since  there  is  less  time  for  reaction 
and  relaxation  to  establish  a  steady  state  before  the  next  incremental 
H2  pressure  change.  The  HFss(v)  increase  and  decay  is  proportional  to 
H2  which  indicates  that  there  is  minimal  F  atom  depletion  and  that 
collisional  relaxation  is  not  important  (see  equation  III-l).  The  decay 
of  the  HFss(v)  begins  when  the  Hj  is  turned  off.  For  t  <  tQn,  the 
"real"  F  is  much  greater  than  Fs*  (less  reaction  than  ocurrs  for  steady 
state)  and  therefore  the  "real"  HF(v)  is  greater  than  HFss(v).  When  the 
H2  is  turned  off,  the  "real"  F  still  decays  (reaction)  and  tries  to 
approach  a  steady  state.  Because  the  H2  is  changing  more  slowly  during 
its  pump  out  stage,  the  F  eventually  reaches  a  steady  state  and  at  this 
time  starts  to  increase.  Likewise,  since  the  F  is  still  reacting  after 
the  H2  is  turned  off,  the  "real"  HF(v)  is  still  being  formed  and 
approaches  its  steady  state  value  at  which  time  it  decays.  Therefore, 
because  the  "real"  CM-PP  system  deviates  from  a  steady  system,  the  HF(v) 
pressure-time  profiles  do  not  decay  until  after  the  H2  is  turned  off. 

Fluxes  are  listed  in  Tables  III-4  thru  III-6  for  low  and  high  flow 
conditions.  For  low  flows,  the  net  flux  shows  a  single  sign  change, 
that  occurs  after  the  H2  is  turned  off,  indicative  of  the  single  maximum 
observed  in  the  HF(v)  pressure-time  profiles.  Depletion  of  HF(v)  is 
predominantly  due  to  flow  from  the  observation  volume  followed  by 
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radiative  relaxation.  Collisional  relaxation  is  not  important. 

Approximately  80%  of  HF(v)  depletion  is  due  to  flow  from  the 
observation  volume/  while  depletion  due  to  radiative  relaxation  is  "25% 
for  HF(v=3)/  "20%  for  HF(v=2),  and  "12%  for  HF(v=l). 

Reaction  is  the  dominant  input  for  all  vibrational  levels. 

Radiative  relaxation  of  HF(v=3)  accounts  for  <10%  of  the  input  to 
HF(v«2).  For  HF(v=l),  radiative  relaxation  from  HF(v=2)  becomes  of 
comparable  importance  as  input  only  after  2  msec. 

The  HF (v)/HF (v=2) -time  plots  (Figure  III-10)  show  linear  behavior  to 
0.5  msec.  One  hundred  percent  of  the  tQn  data  can  be  used  to 
extrapolate  to  zero  time  and  obtain  nascent  distributions. 


2b.  High  Flows  (F/H2  >  0.75),  tQn  »  0.5  msec 

The  time  dependence  of  the  partial  pressures  of  F  and  H2  are  similar 
to  that  for  tQn  =  10  msec  and  tcyc  =  100  msec.  The  F  atoms  are  the 
limiting  reagent  and  are  "80%  depleted  within  0.5  msec.  This  depletion 
is  due  predominately  to  reaction.  The  HF(v)  pressure-time  profiles 
exhibits  a  double  maxima;  the  first  occuring  before  the  H9  is  turned 
off. 

The  HFss(v)-H2  pressure  profiles  show  similar  behavior  as  the 
profiles  in  Figure  III-8  and  this  behavior  is  explained  in  the 
same  way  as  for  the  ton  =  10  msec  case  by  considering  equation  I1I-1. 

The  "real"  F  and  HF(v)  show  larger  deviations  from  steady  state  for  t  < 
tQn  than  was  observed  for  a  t0„  =  10  msec.  However,  because  F  is 
depleted  during  tQ  ,  at  times  >  ton,  the  "real"  F  is  approximately  equal 
to  Fss.  Thus,  the  behavior  of  the  HF(v)  pressure-time  profiles  for  tQn 
=  0.5  msec  and  tQn  =  10  msec  is  essentially  the  same. 

Tables  III-4  thru  III-6  show  the  fluxes  for  high  flow  conditions. 

H2  collisional  relaxation  accounts  for  up  to  70%  of  the  depletion  of 
HF(v)  while  radiative  relaxation  is  less  than  10%. 

Reaction  is  the  dominant  input  to  HF(v=2).  For  HF(v=l),  reaction  is 
dominant  to  0.1  msec  after  which  collisional  relaxation  from  HF(v=2) 
accounts  for  85%  of  the  input.  Radiative  relaxation  from  HF(v=2) 
accounts  for  10%  of  the  input  until  after  2  msec  when  it  increases  to 
"30%. 


The  HF(v)/HF(v=2)-time  plots  are  shown  in  Figure  III-ll.  These 
plots  are  linear  out  to  0.1  msec  (where  relaxation  processes  are 
minimal)  and  when  extrapolated  to  zero  time  give  the  correct  nascent 
distributions.  In  this  case  24%  of  the  ton  data  can  be  used  for  the 
extrapolation;  for  tQn  =  10  msec  only  16%  of  the  tQn  data  could  be  used. 


Comparison  With  CM-PP  Experiment 


Figure  III-12  shows  experimental  HF(v)  population-time  profiles  and 
ratio  plots  for  the  F  +  H2  reaction  using  a  tQn  =  0.5  msec  and  tc  c  =  25 
msec.  These  HF(v)  population-time  profiles  show  broad  maxima  thax  occur 
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at  times  closer  to  when  the  H2  is  turned  off  than  do  the  simulations. 

The  ratio  plots  show  linear  behavior  and  when  extrapolated  to  zero  time 
give  nascent  distributions  in  agreement  with  literature  values. 

Parameters  in  the  model  were  varied  in  order  to  qualitatively 
reproduce  the  experimental  results.  It  was  found  that  by  increasing  the 
flows  of  F  and/or  H2  (increased  [F]  and/or  [ H2 J )  that  the  HF(v)  maxima 
moved  closer  to  the  time  when  the  H2  was  turned  off.  By  increasing  the 
flows  of  reactants,  percent  reaction  is  increased.  The  [F]  and  HF(v) 
reach  their  steady  state  values  at  earlier  times,  but  still  after  the  H2 
is  turned  off.  Increasing  the  k_  rate  constant  for  flow  of  HF(v)  from 
the  observation  volume  also  caused  the  HF(v)  maxima  to  move  closer  to 
the  times  when  the  H2  is  turned  off.  By  increasing  this  rate  constant, 
the  depletion  flux  becomes  dominant  at  an  earlier  time.  As  long  as  the 
IF]  was  in  excess  of  [H2]/  the  ratio  plots  showed  behavior  similar  to 
that  of  Figure  III-10,  otherwise  relaxation  processes  became  important 
and  the  ratio  plots  resembled  those  of  Figure  III-ll. 

The  difference  between  the  experimental  and  simulated  profiles  may 
be  some  combination  of  larger  [F]  and/or  [H2]  and  faster  flow  from  the 
observation  volume.  The  behavior  of  the  experimental  profiles  could  also 
indicate  that  the  HF(v)  reenters  the  observation  zone  after  being 
partially  relaxed  by  either  the  vessel  walls  or  by  gas  phase  collisions. 
As  pointed  out  earlier,  this  path  is  not  included  in  the  simulation. 
However,  if  this  process  occured  it  would  be  expected  that  the  product 
energy  distributions  obtained  would  be  relaxed. 


CM-PP  Conclusions 

There  are  two  limiting  cases  of  experimental  conditions:  low  and 
high  flow  of  H2.  The  experiments  performed  under  low  flow  conditions 
are  characterized  by  a  small  amount  of  F  depletion  and  a  single  maximum 
in  the  HF(v)  pressure-time  profiles.  Collisions  with  H2  are  not  a 
dominant  process  under  these  conditions.  Plots  of  the  pressure  ratios 
of  HF(v)/HF(v=2)  as  a  function  of  time  show  linear  behavior  at  times 
when  relaxation  processes  are  minimal  and  extrapolation  to  zero  time 
gives  the  nascent  product  energy  distributions.  By  decreasing  the 
observation  time,  a  larger  portion  of  the  data  (out  to  the  time  when  the 
H2  is  turned  off)  can  be  used  in  the  extrapolation.  This  is  an 
important  experimental  consideration  because  the  earliest  data  will  have 
the  lowest  intensity  and  thus  the  lowest  signal/noise  ratio.  It  was 
also  shown  that  in  order  to  obtain  nascent  distributions,  the  cycle  time 
must  be  long  enough  to  ensure  that  the  reaction  vessel  is  evacuated 
before  the  next  pulse. 

Double  maxima  are  observed  in  the  HF(v)  pressure-time  profiles  for 
high  flow  conditions.  The  double  maxima  are  due  to  F  atom  depletion  and 
H2  collisional  relaxation.  The  times  at  which  the  ratio  plots  exhibit 
linear  behavior  are  much  shorter  than  for  low  flow  conditions  even  when 
the  observation  time  is  reduced.  This  is  because  H2  collisional 
relaxation  changes  the  relative  composition  of  the  reaction  mixture 
before  the  H2  is  turned  off.  Although  the  nascent  distributions  can  be 
obtained  from  the  simulated  data,  experimentally  it  would  be  more 
difficult  because  of  the  small  signal/noise  ratio  of  these  initial 
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The  "best"  CM-PP  experiment  is  then  one  in  which:  the  observation 
time  is  shorter  than  the  fastest  relaxation  process,  the  cycle  time  is 
long  enough  to  evacuate  the  observation  volume,  a  low  flow  of  H2 
(minimal  F  depletion),  so  that  the  HF(v)  profile  exhibits  a  single 
maximum. 


CM-LP  MODEL  CALCULATIONS 


The  CM-LP  experiments  consist  of  flowing  mixtures  of  SFg  in  Ar,  HR 
(HR  is  a  generic  chemical  entity  attached  to  H)  in  Ar,  Ar,  and  a 
deactivator,  M,  into  a  reaction  vessel.  Fluorine  atoms  are  formed  on  a 
microsecond  time  scale  by  the  multiphoton  dissociation  of  SFg42.  The  F 
atoms  react  with  HR  to  produce  HF(v)  which  can  be  col  1  isional ly 
deactivated,  undergo  radiative  emission,  or  flow  out  of  the  observation 
volume.  The  CM-LP  experiments  are  further  simplified  by  having  [F]  << 
(HR]  (pseudo  first  order  conditions).  For  the  experiments  described  in 
this  report  M  =  HR  although  essentially  any  deactivator  could  be  used. 


The  following  reaction  scheme  is  used  to  model  the  experiments: 

nhv 

SFg  - >  F 

^  i 

F  +  H2  - >  HF  ( i)  +  H 

k j , i (H2) 

HF ( i)  +  H2  —A1 - >  HF( j)  +  H2 

kj  i(HF0) 

HF  (  i)  +  HF  ( 0 )  --- - >  HF  (  j  )  +  HF  (k ) 

ki  jfAr) 

HF  ( i)  +  Ar  -1L - >  HF  ( j )  +  Ar 

Ai  i 

HF  (  i)  - AL - >  HF  ( j )  +  h 


P 

HF  (  i)  - £ - > 

An  equilibrium  (Boltzmann)  rotational  distribution  is  assumed, 
justifiably,  since  these  experiments  are  done  at  total  pressures  of 
approximately  4  torr. 


The  reaction  scheme  can  be  described  by  a  set  of  coupled 
differential  equations  in  which  it  is  assumed  that  i  =  j  +/-  1  (only 
nearest  neighbor  transitions  are  allowed): 

d  [F]/dt  =  -ktot[H2]  [F] 

d[HF(i)]/dt  =  ki[H2](F]  +  kif j(H2) (HF(j) ] [H2] 

-  kjf i(H2)  [HF  (  i ) I  [H2]  -  Ajfi[HF(i)] 

+  A±r  j [HF ( j) ]  -  kj#i(Ar)  [HF(i)] (Ar] 
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+  kif  j  (Ar)  tHF(j)  ]  [Ar] 

-  kjf i(HF0)  [HF (  i) ]  [HF  ( 0 ) ] 

+  kif  j(HF0)  [  HF  ( j )  ]  (HF(O)]  -  kp  [HF  (i)  ] 

These  equations  are  numerically  integrated  using  an  IMSL 
subroutine  (Runge-Kutta-Verner  5th  and  6th  order  method)  to 
obtain  the  time  dependence  of  E  and  HF(v).  A  tolerance  of  0.001 
and  time  interval  of  5.0  x  10“'  sec  was  used  for  the  numerical 
integration.  The  parameters  used  were  the  same  as  for  the  CM-PP 
model  except  for  the  following: 

[ H2 3  =  20-60  mtorr 

[FJq  =  0.1  mtorr 

[Ar]  =  4.0  torr 

k^fj(Ar)  =  30  torr-1  sec-1  (reference  43) 
kp  =  60  sec-1 

Analysis  of  Rate  Constants 

For  the  highest  vibrational  level  populated,  vmax,  the  input  is  due 
solely  to  reaction  and  the  analytical  solution  for  the  time  dependence 
of  HF(vmax)  is  the  sum  of  two  exponentials44: 

HF(v=vmax)  =  Aexp(-kfcott)  +  Bexp(-krf vmaxt)  III-2) 

ktot  =  ktot^H2^ 

kr,  vmax  =  ki,j(H2HH2]  +  kp  +  *i#j  +  (Ar)  [Ar] 

+  kif j(HF0) [HF ( 0 ) ] 

A  =  -B  =  (kvmax  ^h2 HF]  q)/ (kr ,  vmax-ktot  ^ 

The  analytical  solutions  for  the  lower  vibrational  levels  are 
complicated  by  additional  input  due  to  vibrational  relaxation  of  the 
higher  levels.  In  general,  for  each  vibrational  level  below  vmax  an 
additional  exponential  term  appears  in  the  analytical  solution44. 

The  pseudo  first  order  rate  constants  are  obtained  from  a  non-linear 
least  squares  fit  of  the  data  to  these  analytical  solutions.  The  non¬ 
linear  least  squares  fit  was  done  using  the  CURFIT  subroutine  by 
Bevington45  which  minimizes/?=  ^[HFca]-c  (t)  -  HFexP(t)]2/  ^Hp(t)2  with 
respect  to  the  parameters  (pseudo  first  order  rate  constants;.  If  the 
pseudo  first  order  rate  constants  for  relaxation  are  then  plotted  as  a 
function  of  H2  pressure,  the  slope  is  the  H2  collisional  deactivation 
rate  constant,  k^  and  the  intercept  is  the  sum  of  the  remainder  of 

the  rate  constant^ Jf or  depletion  of  HF(v).  Similarily,  a  plot  of  the 
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pseudo  first  order  reaction  rate  constants  versus  H2  pressure  has  a 
slope  equal  to  the  total  reaction  rate  constant,  ktofc. 

Using  the  double  exponential  function,  the  rate  constants  obtained 
for  HF(v=3)  were  in  agreement  with  the  ones  used  in  the  model.  The  rate 
constants  obtained  for  HF(v=2)  and  HF(v=l)  using  the  multiexponential 
function  were  low  by  25%.  For  HF(v=3),  only  4  parameters  must  be  fit  by 
the  non-linear  least  squares,  while  fitting  HF(v=2)  and  HF(v=l)  data 
involves  6  and  8  parameters  respectively. 

An  alternative  method  to  obtain  the  rate  constants  is  to  decouple 
the  HF ( v)  differential  equation  into  components:  relaxation  from  the 
upper  levels  into  v  and  input  to  v  due  to  reaction,  thereby  the  lower 
levels  have  an  analytical  solution  equivalent  to  equation  III-2.  Figure 
III-13  shows  a  two  level  system  (using  the  F  +  H2  reaction  (vmax  =  3) 
as  an  example)  from  which  the  algorithm  for  this  method  is  illustrated. 
For  HF (v=3) ,  input  is  due  only  to  reaction  while  output  is  the  sum  of 
all  the  depletion  processes.  This  can  be  fit  to  the  analytical  solution 
in  equation  III-2.  For  HF(v=2),  input  is  due  to  reaction  and  some 
fraction,  X3,  of  the  population  from  HF(v=3)  that  is  relaxed  into  v=2. 
The  HF(v=2)  population  formed  by  reaction  can  undergo  relaxation 
processes  with  a  pseudo  first  order  rate  constant  (kr  2)  while  the 
HF(v=2)  population  formed  by  relaxation  from  HF(v=3)  can  also  undergo 
relaxation  processes,  k'r  The  terms  kr  2  and  k'  ?  are  equivalent, 
however  the  prime  denotes' the  source.  From  this  model,  a  set  of 


differential  equations  can  be  written: 

dHF ( 3 ) /dt  =  k3[H2][F]  -  kr^  3 [HF (3) ]  III-3) 

dHF(2)/dt  =  {  k  2  [  H2  J  [F  ]  -  krf2[HF(2)]  } 

+  { X 3 k  r  ^  3  [HF  (3)  ]  -  t,2  [HF  (2)  ]  '}  III-4) 

dHF(2)tot/dt  =  dHF (2) rxn/dt  +  dHF(2)'/dt  III-5) 

where  the  {  }  terms  in  equation  III-4  are  replaced  by  the  respective 
dHF(2)/dt  terms  in  equation  III-5.  Integration  of  equation  III-5  for 
the  time  dependence  of  HF(2)  gives: 

HF (2) tot  =  HF { 2 ) rxn  4  HF (2)  *  III-6) 

Equation  III-6  can  be  rearranged  to  give: 

HF (2) rxn  =  HF(2)tot  -  HF  (2) *  III-7) 


Thus  after  subtracting  the  contribution  due  to  HF(v=3)  relaxation  from 
HF(2)t0  ,  the  resulting  HF(2)rxn  has  the  analytical  solution  of  equation 
III-2.  An  iterative  method  is  defined  in  which  a  k'  2  is  estimated, 
dHF(2)'/dt  is  numerically  integrated  to  give  HF(2)','and  HF(2)rxn  is 
determined  from  equation  III-7.  The  time  dependence  of  HF(2)rxn  is 
determined  by  a  non-linear  least  squares  fit  to  equation  III-2  and  a  new 
kr  2  is  obtained.  Convergence  occurs  when  kr  0  is  within  1%  of  the 
estimated  k'r  2«  This  procedure  is  performed  tor  0<X3<1.  The  "best"  fit 
for  a  particular  X3  is  chosen  to  be  the  one  with  the  minimum  chi  square 
statistic.  In  addition,  kfcot  can  be  used  as  a  check  of  self 


consistency,  since  it  is  constrained  to  be  the  same  for  all  vibrational 
levels.  The  iterative  procedure  can  then  be  repeated  for  each  lower 
vibrational  level. 

The  pseudo  first  order  rate  constants  for  reaction  and  relaxation 
obtained  are  then  plotted  as  a  function  of  H2  pressure.  The  slopes  of 
these  plots  give  the  total  reaction  rate  constant  and  deactivation  rate 
constants.  Depending  on  the  relative  magnitudes  of  k^ot  and  k^  ^(M),  the 
pseudo  first  order  rate  constant  obtained  from  the  rise  of  the'HF(v) 
time  profile  can  be  either  due  to  reaction  or  relaxation.  The  same  is 
true  for  the  pseudo  first  order  rate  constant  obtained  from  the  fall  of 
the  HF(v)  time  profile.  Simulations  showed  that  when  k^  ^  (H2 )  <  ktQt, 
the  pseudo  first  order  rate  constant  obtained  from  the  rise  of  the  HF(v) 
time  profile  is  equal  to  kfcot[H23  and  that  from  the  fall  is  due  to 
k^  -jfHoHHo]  and  all  other  deactivation  terms  that  are  independent  of 
H2 •  When  -5(^2)  >  ktot,  the  reverse  is  true.  This  distinction  must 
be  made  in  e'xperiments  where  the  "correct"  process  is  not  known  apriori; 
additional  diagnostic  experiments  must  be  performed. 


CM-LP 


The  F  atoms  decay  exponentially  with  a  pseudo  first  order  rate 
constant,  ktot  [H2 ] -  Figure  III-14  shows  the  HF(v)  pressure-time 
dependence  which  are  observed  in  the  new  CM-LP  experiments.  The  HF(v) 
pressure-time  profiles  show  an  initial  rise  followed  by  a  slower  decay. 
The  decay  is  due  to  collisional  and  radiative  relaxation  and  appears  to 
decrease  for  the  lower  vibrational  levels  because  of  the  input  they  are 
receiving  from  the  upper  levels.  The  HF(v=Q)  shows  a  much  slower  decay 
which  is  due  solely  to  the  k  term.  The  total  emission  intensity,  which 
has  been  observed  in  earl ier  ^experiments,  is  a  composite  of  the  HF(v=l- 
3)  populations  of  Figure  III-14  weighted  by  their  respective  Einstein 
coefficients.  To  be  emphasized  is  that  the  total  intensity-time  profile 
exhibits  a  non-exponential  decay. 

The  time  dependence  for  the  population  ratios  (relative  to  HF(v=2) 
is  shown  in  Figure  III-15.  The  ratio  plots  behave  linearly  out  to  0.3 
msec  and  when  extrapolated  to  zero  time  give  the  nascent  distributions. 

Table  III-7  summarizes  the  rate  constants  obtained  by  analyzing  the 
simulated  HF(v)  popul ations-time  profiles.  The  reaction  rate  constant 
obtained  from  the  simulated  total  intensity  fit  is  approximately  20% 
lower  than  that  input  to  the  model,  while  the  deactivation  rate  constant 
is  approximately  35%  too  low.  Thus  within  the  normal  error  of  these 
experiments,  the  total  intensity  should  not  be  used  to  calculate  these 
rate  constants.  The  v=3  fit  produce  reaction  and  deactivation  rate 
constants  equivalent  to  that  input  to  the  model  as  does  the  fit  to 
HF(v=2)  when  Xv  =  1.0.  The  dependence  of  the  reduced  chi-square 
statistic  on  Xo  is  significant  for  the  simulated  data;  its  ability  to 
differentiate  different  values  of  X2  will  decrease  for  data  that 
includes  random  experimental  error.  For  experimental  data  both  chi- 
square  and  k|-0£  can  be  used  as  a  test  for  the  goodnes  of  fit.  A 
limitation  of  the  iterative  procedure  is  that  when  the  amount  to  be 
subtracted  in  equation  III-6  is  >  70%  of  the  total  amount  present,  the 
algorithm  either  does  not  converge  (2500  iterations  are  performed)  or  it 
produces  incorrect  rate  constants.  This  was  observed  for  HF(v=l) 


RELATIVE  POPULATION 


Table  1 1 1-7 


Rate  Constagts  f or^Simulated  Data 
( x  10°  torr  sec  1 ) 


Vibrational  Level 

k,  . 

tot 

X 

V 

kv-l,v(H2> 

x2(») 

3 

1.00 

— 

0.201 

2 

0.838 

0.0 

0. 149 

7  46  x  101 

2 

0.915 

0.5 

0. 175 

2.22  x  101 

2 

Total  Intensity 

1.01 

0.800 

1.0 

0. 195 

0.130 

3.08  x  10'1 

aBased  on  seven  significant  figures  for 
single  precision  on  the  Prime  computer. 
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simulated  data  where  the  contribution  from  HF(v=2)  relaxation  is  large. 
Another  possible  method  for  analysis  is  to  perform  a  simultaneous  least 
squares  fit  of  the  data  for  HF{v=3),  HF(v=2)  and  HF(v=l)  and  calculate 
all  of  the  rate  constants. 


CM-LP  Conclusions 

CM-LP  experiments  can  provide  nascent  distributions,  total  reaction 
rate  constants,  and  microscopic  deactivation  rate  constants.  The 
analysis  of  these  experiments  are  more  direct  than  the  CM-PP  experiments 
because  of  the  pseudo  first  order  conditions. 

Nascent  distributions  are  obtained  by  extrapolating  time  dependent 
HF(v)  population  ratios  to  zero  time.  Pseudo  first  order  reaction  rate 
constants  and  deactivation  rate  constants  are  obtained  by  performing 
experiments  at  various  H2  pressures  and  fitting  the  time  dependent  HF(v) 
populations  to  the  sum  of  two  exponentials.  For  lower  vibrational 
levels  (<  vmax)  the  contribution  from  relaxation  of  upper  vibrational 
levels  destroys  the  double  exponential  function.  It  was  shown  that  the 
rate  constants  can  be  calculated  by  an  iterative  procedure  when  the 
relaxation  contribution  to  the  input  is  <  70%  of  the  total  population  of 
the  specific  vibrational  level.  The  pseudo  first  order  rate  constants 
are  then  plotted  as  a  function  of  H2  pressure  and  the  slopes  of  these 
plots  give  the  respective  total  reaction  rate  constant  and  microscopic 
deactivation  rate  constant. 


It  was  also  shown  that  fitting  total  intensity-time  profiles  to  the 
sum  of  two  exponentials  produces  rate  constants  that  are  consistently 
low  and  outside  of  the  normal  error  range  for  this  type  of  experiment. 


IV.  EXPERIMENTAL  APPARATUS  and  DATA  PROCESSING 


EXPERIMENTAL  APPARATUS 

Modeling  of  the  first  generation  CM-PP  exper iments14,26a  (see 
section  III)  provided  information  for  the  following  design  improvements: 

1.  Replacement  of  diatomic  reactant  injector  with  a  piezoelectric  valve 
thereby  reducing  the  pulse  length  from  10  msec  to  0.5  msec 

2.  Replacement  of  a  six  channel  signal  averager  (10  sec/channel)  with  a 
2048  channel  signal  averager  (0.5  microseconds/channel) 

3.  The  White  cell  is  replaced  with  a  focusing  lens  and  rear  mirror 

4.  The  1500  cc  reaction  vessel  was  replaced  with  a  10  cc  reaction  vessel 

Figure  IV-1  displays  the  modified  CM-PP  apparatus  while  Figure  IV-2 
shows  the  CM-LP  apparatus. 

A  CM  experiment  provides  both  spectral  and  temporal  resolution. 
Infrared  emission  from  the  products  formed  by  reaction  of  an  atomic 
species  and  hydrogen  donor  molecule  is  focused  by  a  lens  onto  the  slit 
of  a  monochromator  and  dispersed  by  a  grating  (30C  lines/mm);  total 
intensity  is  measured  using  a  mirror  in  place  of  the  grating.  A 
microcomputer  controls  the  stepping  motor  that  moves  the  grating  through 
the  wavelength  region  of  interest  (2. 5-5.0  microns).  At  each 
wavelength,  an  emission  intensity-time  profile  is  measured.  The 
infrared  emission  is  detected  by  a  liquid  nitrogen  cooled  InSb  detector, 
amplified,  collected  by  either  a  transient  recorder  or  signal  averager 
and  stored  on  a  floppy  disk  from  which  it  is  transferred  to  a  mainframe 
computing  system  for  further  data  analysis.  The  individual  parts  of  the 
CM-PP  and  CM-LP  experimental  apparatus  will  be  discussed  separately  in 
the  following  sections. 


CM-PP  Experimental  Apparatus 


Flow  System 

Tank  gases  were  used  from  cylinders  with  regulators  attached  while 
CH2CI2  was  stored  in  a  130  cc  glass  bulb  with  stopcock  and  was  degassed 
by  the  f reeze-pump-thaw  method  prior  to  use.  Copper  lines  attached  to 
the  cylinder  regulators  or  to  the  CH2Cl2  glass  bulb  allow  the  gases  to 
flow  into  a  glass  high  vacuum  system  which  leads  to  the  reaction  vessel. 
Separate  flow  lines  with  Helicoid  pressure  gauges  and  greaseless 
stopcocks  were  used  for  the  atomic  precursor  and  hydrogen  (deuterium) 
donor  reagent.  Fluorine  and  hydrogen  bromide  were  handled  in  passivated 
copper  lines  attached  directly  to  the  reaction  vessel.  These  gases  were 
stored  and  handled  in  a  well  ventilated  hood.  All  reagent  (except 
fluorine)  flows  were  controlled  by  MRS  electronic  flowmeters  (model  258) 
with  MRS  control  valves  (type  248A)  and  an  MRS  control  system  (series 
260).  Fluorine  flows  were  controlled  by  a  needle  valve  and  monitored 
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Figure  I V- 1 .  Diagram  of  CM-PP  Experimental  Apparatus 


with  the  Helicoid  pressure  gauge. 

A  Cooke  cold  cathode  vacuum  gauge  was  used  to  measure  the  pressure 
near  the  pumpout  of  the  reaction  vessel.  A  Cooke  thermocouple  guage 
measured  pressure  along  the  discharge  tube.  Mechanical  Welch  Duo  Seal 
vacuum  pumps  were  used  to  evacuate  the  flow  lines.  The  reaction  vessel, 
a  4"  aluminum  cube  with  1"  holes  drilled  thru  the  center  of  each  face, 
was  pumped  out  by  a  Welch  Duo  Seal  pump  in  series  with  a  6"  liquid 
nitrogen  trapped  oil  diffusion  pump  {NRC  VHS-6).  A  6"  slide  valve 
separates  the  oil  diffusion  pump  from  the  cold  nitrogen  trap  attached  to 
the  reaction  vessel.  Effluent  from  the  reaction  vessel  was  vented  to 
the  hood. 

Reagent  Injection 

A  Lasertechnics  Model  LPV  piezoelectric  valve  was  used  to  inject  the 
hydrogen  donor  reactant  into  the  reaction  vessel.  A  Lasertechnics  Model 
203  LPV  valve  driver  provides  electrical  pulses  to  control  the 
piezoelectric  valve.  For  the  reported  experiments,  the  valve  was 
operated  with  a  pulse  width  of  0.5  msec  and  a  cycle  time  between  25  and 
30  msec  (pulse  rate  of  33-40  Hz).  A  nozzle  with  a  1.0  mm  cylindrical 
orifice  was  used  to  minimize  cluster  formation.  The  valve  driver 
triggers  the  signal  averager  to  start  data  collection.  The  pulse  width 
and  delay  between  pulses  is  determined  and  checked  during  the  experiment 
by  connecting  the  trigger  output  from  the  valve  driver  to  a  Tektronix 
533A  oscilloscope.  The  cycle  time  is  varied,  depending  on  the 
experiment,  to  ensure  that  the  vessel  is  completely  pumped  out  (no 
emission  observed)  before  the  next  pulse. 


Atom  Generator 

The  microwave  discharge  of  CF^  gas  generated  F  atoms  while  hydrogen 
or  deuterium  gas  was  used  to  produce  H  and  D  atoms  respectively.  The 
microwave  power  was  supplied  by  a  115  watt  2450  MHz  Raytheon  model  PGM 
10X2  microwave  generator  coupled  to  a  Raytheon  Model  E14C  Evenson  type 
microwave  cavity.  The  atomic  precursor  flows  through  a  12  mm  o.d. 

Vycor  discharge  tube  located  along  the  axis  of  the  cavity;  the  exit 
orifice  is  1  mm  in  diameter.  The  discharge  tube  was  inserted  into  the 
reaction  vessel  (approximately  0.5  cm  from  the  center  of  the  vessel) 
using  a  Cajon  Ultra  torr  O-ring  seal  and  fitting.  Air  flowing  through 
the  cavity  efficiently  cools  the  discharge  region.  When  CF^  was  the 
atomic  precursor,  a  clean  Vycor  discharge  tube  was  used  which  needed  to 
be  replaced  after  approximately  8  hours  of  CF*  flow  due  to  etching  by  F 
atoms.  The  inside  surface  of  the  discharge  tube,  when  used  with 
hydrogen  or  deuterium,  was  treated^®  prior  to  use  in  order  to  decrease 
atom  recombination  on  the  walls.  As  in  reference  46,  the  tube  was 
filled  with  a  10  M  NaOH  solution  for  approximately  18  hours  followed  by 
several  washings  with  distilled  water.  The  tube  was  then  filled  with  a 
10  M  HNO3  solution  for  approximately  18  hours,  washed  several  times  with 
distilled  water  and  dried. 

Detection  System 


Infrared  radiation  from  the  reaction  products  is  focused  by  a  CaF2 


focusing  lens  (diameter  of  2"  and  focal  length  of  6")  onto  the  slit  of  a 
1/4  meter  near  Littrow  monochromator  (Perkin  Elmer  model  99G)  which  is 


coupled  to  a  liquid  nitrogen  cooled  InSb  detector  (Infrared 
Laboratories)  by  a  pair  of  ofi  axis  ellipsoidal  mirrors.  A  NEP 
equivalent  ^>ower)  of  1  x  10“13  watt/dizl1/”  at  2  microns  or  4  x 

wat-fc/fHlzl  -*■'  ^  5  mirrnns  use  moa  rii  >  4  '  .  The  o/->4- n  r  recnr 
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equivalent  ^>ower)  of  1  x  lO-1-3  watt/(HzV L/  “  at  2  microns  or  4  x  10-14 
watt/(Hz) r  at  5  microns  was  measured39'47.  The  detector  response  was 
<  1  microsecond  when  using  a  10K  feedback  resistor.  The  wavelength 
sensitivity  of  the  monochromator-detector  system  was  calibrated  with  a 
blackbody  radiation  source  (Infrared  Industries  model  463)  and 


temperature  controller  (model  101B). 


In  the  previous  CM  experiments,  a  2. 5-5. 5  micron  spectral  bandpass 
filter  was  positioned  between  the  InSb  detector  and  an  external  KRS-5 
vacuum  window  attached  to  the  dewar  housing.  Both  the  filter  and 
detector  were  cooled  to  liquid  nitrogen  temperatures,  thereby  reducing 
the  background  radiation  and  thermal  noise  viewed  by  the  detector. 
Blackbody  measurements  using  this  configuration  showed  the  presence  of 
an  absorption  between  2.9  and  3.3  microns  that  became  enhanced  (larger 
and  wider)  over  a  2  hour  time  period.  In  this  work,  a  room  temperature 
filter  (2. 5-5. 5  microns)  positioned  between  the  exit  slit  of  the 
monochromator  and  the  KRS-5  vacuum  window  replaced  the  cold  filter. 
Blackbody  measurements  indicated  that  the  absorption  was  eliminated  but 
with  a  reduction  in  signal/noise.  Blackbody  measurements  also  showed 
absorptions  between  2.6  and  2.9  microns  and  between  4.2  and  4.4  microns 
which  were  due  to  H20  (2.7  microns48)  and  C02  (4.3  microns48) 
respectively.  These  absorptions  were  eliminated  when  the  optical  path 
was  enclosed  by  a  cylinder  through  which  dry  air  or  nitrogen  (for  the  D 
+  F,  experiments)  was  purged  for  approximately  8  hours  prior  to 
performing  the  experiment. 

The  signal  from  the  detector  was  ac  coupled  to  an  EG&G  model  113 
preamplifier.  For  the  CM-PP  experiments,  the  high  frequency  rolloff  was 
set  at  30KHz  which  was  the  optimum  setting  to  eliminate  high  frequency 
noise  without  distortion  of  the  experimental  signal. 

An  EG&G  model  4203  signal  averager  was  used  for  data  collection. 

When  the  signal  averager  was  triggered  by  the  pulsed  valve  driver,  it 
started  data  aquisition.  The  signal  averager  stores  2048  data  points 
with  a  time  resolution  of  5  microseconds.  The  experimental  signal  is 
transferred  via  a  customized  parallel  interface  to  a  Z80A  based 
microcomputer  for  further  data  analysis.  The  microcomputer  stores  128 
data  points  on  a  floppy  disk,  each  of  which  is  the  average  of  two  of  the 
data  points  from  the  signal  averager.  Due  to  the  inherent  mechanical 
delay  between  the  driver  and  valve,  the  first  10  data  points  stored 
represent  the  baseline  and  are  necessary  to  determine  the  baseline  for 
the  data  analysis. 


A  typical  CM  experiment  consists  of  stepping  through  the  wavelength 
region  of  interest,  the  number  of  steps  being  between  200  and  350 
depending  on  the  experiment,  with  step  sizes  ranging  from  3  to  5  cm~^. 
The  spectral  resolution  is  between  11  and  22  cm-r  and  is  determined  by 
the  monochromator  slit  width  used  (1.0  -  2.0  mm)  in  order  to  obtain 
sufficient  measurable  intensity.  Total  HF(v)  intensity  is  measured 
before  and  after  each  experiment  to  ensure  that  flows  remained  constant. 
The  microcomputer  controls  the  stepping  motor  that  rotates  the  grating 


of  the  monochromator  and  also  controls  transfer  of  data  from  the  signal 
averager. 


CM-LP  Experimental  Apparatus 


Flow  System  And  Reagent  Mixtures 

Mixtures  of  atomic  precursor  (SFg)  or  HR  reactant  in  Ar  were  made 
and  stored  in  12  liter  glass  bulbs  with  stopcocks  and  a  MRS  capacitance 
manometer  (model  222,  1000  torr)  was  used  to  measure  pressures. 

Mixtures  were  made  by  filling  the  12  liter  storage  vessel  first  with  70- 
100  torr  of  HR  or  SFg.  These  gases  were  used  directly  from  their 
cylinders  with  attached  regulators.  Methylene  chloride  was  used  from  a 
130  cc  glass  bulb  with  a  stopcock  and  was  purified  by  the  freeze-pump- 
thaw  method  prior  to  use.  The  12  liter  storage  vessel  was  then  filled 
with  Ar  to  a  total  pressure  of  approximately  700  torr.  Argon  was  used 
from  its  cylinder  with  attached  regulator  and  flowed  through  a  liquid 
nitrogen  cooled  molecular  sieve  (Union  Carbide  type  AW-500  pore  diameter 
5  A)  trap  to  remove  any  H2O.  Typically  10%  mixtures  of  HR  in  Ar  and  20% 
mixtures  of  SFg  in  Ar  were  used  approximately  24  hours  after  their 
preparation.  The  12  liter  glass  bulbs  were  connected  to  the  glass  high 
vacuum  system  by  O-ring  seals,  #9  glass  joints  and  Cajon  stainless  steel 
flexible  tubing.  Argon  was  used  from  a  cylinder  with  a  gas  pressure 
regulator  attached  to  copper  lines  leading  to  the  glass  high  vacuum 
system.  The  mixture  flows  were  measured  with  Matheson  #610A  flowmeters 
and  controlled  with  needle  valves.  Argon  flows  were  measured  with  a 
Fischer  &  Porter  flowmeter  and  controlled  with  a  needle  valve.  The  two 
mixtures  and  Ar  flow  into  the  reaction  vessel  with  the  Ar  first  flowing 


through  a  liquid  nitrogen  cooled  molecular  sieve  trap  to  remove  any  H2O. 
Reagent  flows  ensured  that  the  reaction  vessel  was  completely  flushed 
twice  between  laser  pulses  and  that  pseudo-first  order  experimental 
conditions  (F  <<  H2)  are  maintained.  Welch  Duo  Seal  mechanical  pumps 
were  used  to  pump  out  the  system. 


Reaction  Vessel 

The  reaction  vessel  for  these  experiments  was  a  cross  shaped  glass 
tube.  Sodium  chloride  windows  (1"  in  diameter  and  0.25"  thick)  were 
attached  to  each  of  the  four  ends  of  the  tube  via  #25  glass  joints  and 
clamps.  The  vessel  was  attached  to  the  high  vacuum  system  and  pumping 
system  via  ball  and  socket  joints.  The  reagent  mixtures  entered  through 
an  inlet  attached  to  one  arm  of  the  vessel  and  was  pumped  through  an 
outlet  attached  to  the  opposite  arm.  A  COo  laser  was  focused  on  the 
center  of  the  vessel  by  a  BaF2  lens  (focal  length  of  4"  and  1.5"  in 
diameter).  Ninety  degrees  from  the  laser  beam  an  aluminum  coated  mirror 
was  placed  and  opposite  the  mirror  was  a  CaF2  focusing  lens  (diameter  of 
2"  and  focal  length  of  6")  and  detection  system.  The  vessel  and  optical 
path  was  enclosed  in  an  aluminum  box  through  which  dry  air  was  flowed  to 
eliminate  absorption  of  the  2.7  micron  HF  emission  by  H2O  in  the  air.  A 
MRS  capacitance  manometer  (model  222,  10  torr)  was  attached  to  the  top 
of  the  vessel  to  measure  the  pressure. 
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Atom  Generator 


A  pulsed  (1  Hz)  CO?  Tachisto  TEA  laser  (model  400XR)  at  10.59 
microns  (001-100  P{20)  “transition)  was  used  to  generate  F  atoms  by  the 
multiphoton  dissociation  of  SFg.  A  Gentec  joulemeter  (model  ED500)  was 
used  to  measure  the  pulse  energy;  the  output  is  observed  on  a  Tektronix 
533A  oscilloscope.  Typical  pulse  energy  was  0.5  joules.  Within  5%,  the 
pulse  energy  and  profile  were  reproducible  from  pulse  to  pulse. 


Detection  System 

The  detection  system  for  the  CM-LP  experiments  is  the  same  as 
described  for  the  CM-PP  experiments  except  for  the  following  changes. 

The  high  frequency  rolloff  and  gain  of  the  preamplifier  were  100  kHz  and 
10  Hz  respectively.  The  EG&G  4203  signal  averager  was  replaced  with  a 
Biomation  model  805  waveform  recorder  with  time  resolution  of  0.5 
microsecond  and  the  capability  of  pretriggering.  The  trigger  pulse  to 
the  Biomation  was  provided  by  a  photodiode  directed  at  the  laser  cavity. 
Pretriggering  allows  the  collection  of  baseline  points  prior  to 
collection  of  the  signal;  baseline  correction  of  the  experimental 
signal  is  necessary  for  correct  data  analysis.  The  signal  (2048  data 
points)  is  transferred  to  the  microcomputer  via  a  custom  parallel 
interface.  Every  16  data  points  from  the  transient  recorder  are 
averaged  to  give  1  data  point;  thus  a  total  of  128  data  points  are 
stored  temporarily  by  the  microcomputer.  Five  hundred  of  these  data  sets 
are  coadded,  providing  a  factor  of  (500) improvement  in  signal/noise, 
and  finally  stored  on  a  floppy  disk. 

Total  HF (v)  emission  intensity  was  recorded  before  and  after  each 
experiment.  It  was  found  that  the  total  intensity  decreased  by  10% 
during  an  experiment  because  of  the  finite  volume  of  gas  mixtures  used. 

A  linear  correction  based  on  these  total  intensity  measurements  was  made 
to  the  experimental  data.  Laser  noise  affected  only  the  first  data 
point  of  the  spectrally  resolved  emission  intensity-time  profile  and  it 
was  excluded  from  the  data  analysis.  Laser  noise  was  insignificant  in 
the  total  intensity-time  profiles  due  to  the  larger  measured  intensity 
compared  to  the  spectrally  resolved  intensity. 


DATA  ANALYSIS 

In  CM  experiments,  the  intensity  of  individual  vibration-rotation 
states  is  measured  as  a  function  of  both  wavelength  and  time.  Typical 
CM-PP  and  CM-LP  experimental  intensity-frequency-time  surfaces  are  shown 
in  Figures  III-3  and  III-4.  Generally,  individual  vibrational- 
rotational  transitions  are  not  completely  resolved.  Therefore,  in  order 
to  obtain  useful  information  (product  energy  distributions,  rate 
constants)  from  these  spectra,  vibrational-rotational  population 
distributions  must  be  extracted.  The  following  equation1*9  relates  the 
observed  intensity  to  the  populations  of  the  individual  vibration- 
rotation  levels: 


IV-1) 


1(A)  =  Av,jSv,j<»TO0PV(J 

Av  j  =  Einstein  coefficient^®  for  emission  from  vifcurational  level  v 
and  rotational  level  J  in  units  of  sec-1 

Sv  j ( A)  =  Instrumental  response  function  for  a  particular  v,J 
'  transition 

T ( A)  =  Detector  sensitivity 

Pv^j  =  Populations  of  vibration-rotation  levels  v,J  respectively 

I  ()0  -  experimental  intensity  of  a  given  wavelength  in  units  of 
photons/second 

The  instrumental  response  function  (monochromator  slit  function)  is 
represented  as  a  triangle  function: 

sVfJ(A)  -  i-(  A  -  A(v,j»/a 

for  (A-  ( v#  J) )  <  and 

sv, J(A)  -  0 
f  or  (  \  -  X  (v,  J) )  > 

=  monochromator  resolution 

Using  equation  IV-1,  an  iterative  linear  least  squares  is  performed  for 
each  experimentally  measured  time.  In  the  least  squares  method, 

=  <£(KA)calc-  K\)exp)2/^Qj 

is  minimized  with  respect  to  the  vibrational-rotational  populations 
(parameters  of  the  fit).  If  negative  populations  were  obtained  in  the 
fitting  process,  these  levels  were  subsequently  removed  from  the 
analysis  and  another  iteration  performed. 

Since  the  errors  in  the  emission  intensities  were  determined  to  be 
independent  of  wavelength  and  time,  equal  weighting  was  used  in  the 
least  squares  technique.  From  the  vibrational-rotational  populations 
obtained,  a  calculated  spectrum  was  constructed;  the  difference  between 
the  observed  and  calculated  spectrum  was  used  as  an  estimate  of  the 
final  error  in  the  vibration-rotation  populations.  The  rotational 
populations  obtained  from  the  least  squares  fit  exhibited  similar 
distributions  (Boltzmann)  for  each  vibrational  level  and  were  used  as  a 
useful  indicator  that  peak  assignments  were  correct.  The  total 
population  of  a  given  vibrational  level  was  obtained  by  summing  the 
rotational  populations  of  that  level. 

Since  the  rotational  distributions  were  found  to  be  Boltzmann,  an 
alternative  method  of  obtaining  total  vibrational  populations  is  to 
perform  a  nonlinear  least  squares  fit  of  the  data.  The  function  used 
for  this  fit  assumes  a  Boltzmann  rotational  distribution: 


IV-2) 


I(/b  =2"pvav  Zsv  j (  a ) T (A  )  (2J+l)exp(-Ej/kT) 

/■  ~ 

The  nonlinear  least  squares  analysis  was  done  with  the  subroutine  CURFIT 
by  Bevington45  in  which 

-%  =£(KA)calc  -  lWe,P|2/<^) 

is  minimized  with  respect  to  the  vibrational  populations  and  the 
temperaturef  T  (the  parameters  of  the  fit). 

To  obtain  nascent  product  energy  distributions,  ratios  of 
experimental  vibrational  populations  as  a  function  of  time  were 
extrapolated  to  zero  time.  A  weighted  linear  least  squares  fit  was  used 
for  the  extrapolation  where  the  weights  used  were  obtained  from  the 
intensity-population  conversion.  Only  data  out  to  the  maximum 
population  (peak)  was  used  in  the  extrapolation  and  the  zero  of  time  was 
determined  from  the  initial  rise  of  the  total  intensity-time  plots. 

The  analysis  of  HF(v)  population-time  profiles  for  extracting 
reaction  and  deactivation  rate  constants  has  been  described  in  section 
III.  The  data  points  at  long  times  (close  to  baseline  noise)  were 
excluded  from  the  nonlinear  least  squares  fitting  process.  This  was 
done  because  in  the  intensity  to  population  conversion,  the  procedure  of 
removing  negative  populations  causes  these  populations  to  have  only 
positive  contributions  (no  negative  noise). 
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V.  CM  EXPERIMENTAL  RESULTS 


RESULTS  and  DISCUSSION 

Computer  simulations  provided  the  optimum  experimental  conditions 
needed  to  obtain  nascent  product  vibrational  energy  distributions  and 
rate  constants. 

The  results  for  each  chemical  system  are  described  separately 
below.  Typical  pressures  in  the  CM-PP  experiments  measured  at  the  pump 
out  end  of  the  reaction  vessel  were  6  x  10“5  -  1.4  x  10“4  torr.  These 
experiments  were  done  with  a  tQn  =  0.5  msec  and  cycle  times  of  25-30 
msec.  The  total  pressure  measured  in  the  reaction  vessel  for  the  CM-LP 
experiments  was  approximately  4  torr.  The  following  partial  pressures 
were  used  in  the  CM-LP  experiments:  SFg  =  60  mtorr,  HR  =  10-60  mtorr,  Ar 
*  3.5  torr.  Errors  are  reported  at  the  95%  confidence  level. 


F  +  H2 

This  system  has  been  well  studied®'^1**, 26a, 51-54  and  was  used  in 
this  work  as  benchmark  verification  of  both  the  CM-PP  and  CM-LP 
experimental  techniques.  The  total  energy  available  to  the  HF  products 
is  35  kcal/mole5  .  The  third  vibrational  level  lies  32.5  kcal/mole4 
above  the  zero  point  energy  of  the  products  and  is  the  maximum  level 
that  can  be  populated.  In  the  CM-PP  experiments,  CF4  flow  rates  ranged 
from  2.6-46  mole/sec  and  H2  flows  ranged  from  2.0-46  mole/sec.  In  the 
CM-LP  experiments,  the  pressures  used  were:  SFg  =  60  mtorr,  H2  =  20-60 
mtorr,  Ar  =  3.5  torr. 

Typical  HF (v)  total  population-time  profiles  and  HF(v)  ratio-time 
profiles  obtained  from  the  CM-PP  and  CM-LP  experiments  are  shown  in 
Figures  V-l  and  V-2.  The  CM-PP  HF(v)  total  population-time  profiles 
show  a  single  maximum  indicating  the  absence  of  collisional  relaxation 
and  minimal  fluorine  atom  depletion.  When  comparable  CF4  and  H2  flows 
(i.e.  46  moles/sec  for  both  reagents)  were  used,  the  total  population¬ 
time  profiles  peaked  before  the  H2  was  turned  off  and  showed  the  start 
of  a  second  peak  after  the  H2  was  turned  off  indicating  high  flow 
conditions.  These  experiments  were  not  included  in  the  determination  of 
the  reported  initial  distributions.  The  HF(v)  ratio  plots  from  both 
experimental  methods  show  linear  behavior  and  when  extrapolated  to  zero 
time  give  nascent  product  energy  distributions  in  agreement  with 
previously  reported  distributions. 

Reaction  and  deactivation  rate  constants  (Table  V-l)  were  obtained 
from  the  CM-LP  experimental  data.  Figures  V-3  thru  V-7  show  typical 
nonlinear  least  squares  fit  of  the  total  intensity,  and  of  the  HF(v=3) 
and  HF(v=2)  time  profiles  to  the  sum  of  two  exponentials.  Also,  shown 
in  these  figures  are  the  residuals  (HFcalc  -  HFexP)  of  these  fits. 

The  total  reaction  rate  constant  obtained  from  the  analysis  of  the 
total  intensity-time  profile,  0.96  +/-  0.03  x  10&  torr-1  sec-1,  is  in 
agreement  with  other  published  work  in  which  total  intensity 
measurements  are  used  (see  Table  V-l).  From  the  HF(v=3)  population-time 
profiles,  a  reaction  rate  constant  of  1.4  +/-  0.3  x  106  torr-1  sec-1  is 


Reaction^  and  Deactivation  Rate  Constants 
for  the  Reaction  F  +  Hg  - >  HF(v)  +  H 
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Figure  V-5.  Least  Squares  Fit  of  Experimental  HF(v=2) 
Population  Versus  Time,  (X„  =  0.0),  F  +  H0 
Reaction  J  2 


Figure  V-6.  Least  Squares  Fit  of  Experimental  HF(v=2) 
Population  Versus  Time,  (X^  =  0.5),  F  +  H, 
Reaction 


obtained.  As  has  been  shown  by  modeling  in  this  work,  analysis  of  the 
total  intensity-time  profiles  give  consistently  small  values  of  the  rate 
constants  when  forced  to  fit  the  sum  of  two  exponentials,  therefore  the 
result  from  the  fit  of  HF(v=3)  is  preferred. 

The  reaction  rate  constant  obtained  from  the  HF(v=2)  population¬ 
time  profiles  is  1.4  +/-  0.1  x  106  torr-1  sec-1  and  is  in  agreement  with 
that  obtained  from  HF(v=3)  as  expected.  The  quoted  error  is  propagated 
from  the  error  in  the  populations.  Since  error  can  also  be  introduced 
by  the  iterative  method  used  to  obtain  this  rate  constant,  a  more 
realistic  estimate  of  the  error  would  be  a  factor  of  two  larger. 
Likewise,  the  errors  reported  for  the  rate  constants  obtained  from 
analysis  of  the  total  intensity  are  propagated  from  errors  in  the 
emission  intensity  and  do  not  include  the  potential  error  from  using  the 
wrong  function.  Plots  of  pseudo  first  order  reaction  rate  constants 
versus  H2  pressure  for  total  intensity,  HF(v=3)  and  HF(v=2)  gave 
intercepts  of  zero  indicating  that  other  sources  of  F  atom  depletion 
were  absent. 

The  rate  constant  far  deactivation  of  HF(v=3)  by  H2  found  in  this 
work  is  1.5  +/-  0.6  x  105  torr-1  sec"1  which  is  considerably  larger  than 
previously  reported  measurements  (Table  V-l).  A  larger  rate  constant 
suggests  the  presence  of  deactivators  other  than  H2.  In  initial 
diagnostic  experiments  of  this  system  plots  of  pseudo  first  order 
deactivation  rate  constants  versus  H2  pressure  gave  varying  slopes  and 
intercepts,  the  largest  of  which  were  4  x  105  torr-1  sec-1  and  7000  sec- 
1  respectively.  The  intercept  of  these  plots  is  the  sum  of  all 
relaxation  processes  that  are  independent  of  H2  pressure  and  include 
radiative  relaxation,  deactivation  by  Ar  and  HF(v=0),  and  flow  out  of 
the  observation  volume.  When  Ar  was  passed  through  a  liquid  nitrogen 
cooled  molecular  sieve  trap,  the  rate  ccmstant  and  intercept  was  reduced 
to  1.5  x  10^  torr  1  sec-1  and  1200  sec-1  respectively.  These 
diagnostic  experiments  showed  that  H20  was  present  in  the  Argon  and 
although  this  source  of  H20  was  eliminated,  it  is  possible  that  other 
sources  (flow  lines,  vacuum  system,  and  reaction  vessel)  still  remain 
and  are  responsible  for  the  efficient  deactivation. 

F  +  D2 

The  energy  available  to  the  DF  product  is  35 J'  kcal/mole6  allowing 
up  to  the  4th  vibrational  level  (31.69  kcal/mole55  above  zero  point 
energy)  to  be  populated.  This  reaction  is  useful  for  observing  isotope 
effects  and  also  provides  an  additional  verification  of  the  CM-PP 
experimental  method.  The  reaction  forming  DF  has  been  reported  to  be  a 
factor  of  2  slower  than  the  reaction  forming  HF55  and  the  Einstein 
coefficients4  for  DF  emission  are  smaller  (x 4 )  than  for  HF  emission. 
Thus,  higher  reagent  flows  than  were  used  in  the  F  +  Ho  system  were 
necessary  in  order  to  study  this  reaction.  The  flow  of  CF4  was  3.7 
mole/sec  and  D2  flows  were  3.7-10.0  mole/sec. 

The  DF(v)  total  population-time  plots  and  ratio-time  plots  obtained 
from  CM-PP  experiments  are  shown  in  Figure  V-8.  The  nascent 
distribution  (Table  V-2)  obtained  is  in  agreement  with  previous  results, 
peaking  at  the  third  vibrational  level.  Rotational  population  plots  for 
each  vibrational  level  are  Boltzmann  with  T  =  300  K.  The  isotopic 
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Relative  Vibration  Distributions 
for  the  Reaction  F  +  Dg  - >  DF(v)  +  D 
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2 
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0.57 

+/-  0.16 

3 

1.00 

1.00 
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aFrom  reference  7,  value  in  parenthesis  estimated. 

^From  reference  6. 
cFrom  reference  26a. 

^This  work. 

eAverage  fraction  of  the  total  energy  that  goes  to  vibration 
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invariance  of  the  F  +  H2  and  F  +  D2  reaction  is  indicated  by  the  similar 
average  fraction  of  energy  channeled  into  product  vibration,  <fv>,  which 
from  the  CM-PP  results  was  determined  to  be  0.67  and  0.63  for  HF  and  DF 
respectively.  This  efficient  conversion  of  energy  into  product 
vibration  is  due  to  mixed  energy  release  on  a  repulsive  surface56. 


H  +  F2 

The  total  ex^rgy  available  to  the  HF  product  from  this  reaction  is 
101.9  kcal/mole55  allowing  up  to  vibrational  level  11  (99.1  kcal/mole4 
above* exp  point  energy)  to  be  populated.  In  the  previous  CM- 
ppl4,18,26a  msec  F2  "on"  time  and  100  msec  cycle  time)  experiments, 
an  apparent  F2  flow  dependence  was  observed.  When  low  flows  of  F2  (18 
moles/sec)  were  usedj  £he  distribution  peaked  at  v=6  in  agreement  with 
other  studies11'12'57'58  of  this  reaction.  However,  a  population  "hole” 
at  v=4  and  large  populations  in  the  low  vibrational  levels  was  also 
reported  which  had  not  been  observed  before.  For  high  F,  flows  (44 
moles/sec),  a  bimodal  distribution  peaking  at  v=2  and  v=6  was  reported. 
Vibrational  levels  4  and  5  were  highly  depleted.  CM-PP  computer 
simulations  show  that  for  high  flow  conditions,  double  maxima  in  the 
HF(v)  population-time  profiles  should  be  observed  and  that  collisional 
relaxation  is  important.  The  double  maxima  were  probably  masked  in  the 
previous  CM  experiments  since  a  6  channel  signal  averager  with  a  time 
resolution  of  5  msec  was  used.  Collisional  relaxation  would  account  for 
the  large  populations  observed  in  these  experiments  for  the  low 
vibrational  levels.  The  apparent  depletion  of  vibrational  levels  4  and 
5  was  found  to  be  due  to  a  cooled  2. 5-5. 5  micron  spectral  bandpass 
filter  positioned  between  the  InSb  detector  and  a  KRS-5  vacuum  window 
attached  to  the  detector  housing  as  discussed  in  section  IV  of  this 
report.  To  conclude  whether  the  previous  CM  results  were  influenced  by 
the  above  mentioned  factors,  the  H  +  F2  reaction  was  restudied. 

The  results  reported  here  are  from  a  CM-PP  experiment  inwhich  a  H2 
flow  of  3.7  moles/sec  and  a  F2  flow  of  14  moles/sec  was  used.  The 
vibration-rotation  populations  were  obtained  by  using  equation  IV-2  in 
which  a  Boltzmann  rotational  distribution  is  assumed.  The  rotational 
temperature  was  a  variable  parameter  in  the  fitting  process  and  was 
approximately  400  K.  The  assumption  of  a  Boltzmann  rotational 
distribution  is  justified  because  the  H  +  F2  reaction  rate  is  17  times 
slower58  than  the  F  +  H2  reaction  rate  and  tne  reagent  flows  were  higher 
than  used  in  the  F  +  H2  experiment  where  Boltzmann  rotational 
distributions  were  observed. 

The  HF(v)  total  population-time  profiles  and  HF(v)  ratio-time  plots 
are  shown  in  Figures  V-9  and  V-10.  To  avoid  congestion  in  these  plots, 
every  other  vibrational  level  is  plotted.  With  the  exception  of  the 
earlier  reported  CM  results,  the  nascent  product  energy  distribution  is 
in  agreement  with  previously  determined  distributions  (Table  V-3).  For 
the  reasons  given  above,  the  new  distribution  obtained  from  the  CM-PP 
experiments  is  preferred  over  the  previous  CM  results. 

The  <fv>  for  the  H  t-F2  reaction  is  less  than  for  the  F  +  H2 
reaction  while  the  <f^.>11  is  larger.  This  behavior  is  described  by  the 
light  atom  anomaly56  in  which  the  H  atom  comes  to  within  normal  bonding 
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distance  of  the  F2  before  the  repulsive  energy  is  released.  Thus,  a 
larger  fraction  of  the  available  energy  is  channeled  into  relative 
translation  of  the  separating  products. 


D  +  F2 

The  available  energy  for  the  DF  product  is  103.2  kcal/mole55  so  that 
vibrational  levels  up  to  15  (103  kcal/mole4  above  the  zero  point  energy) 
could  be  populated.  Three  previously  reported  experimental  studies  of 
the  D  +  ?2  reaction  are  not  in  agreement  as  to  the  most  populated 
vibrational  level.  The  initial  product  vibrational  energy  distribution 
from  Measured  Relaxation  experiments13  peaks  at  v=10,  fr cm  East  Flow 
experiments59  at  v=9,  while  the  previous  CM  exper iments14,18,26a  gave  a 
narrow  distribution  peaking  at  v=ll  with  no  population  observed  in  v=l- 
9.  The  previously  reported  CM  experiments  were  done  with  high  flows  of 
F 2  (42.6  moles/sec)  and  a  10  msec  observation  time;  conditions  which 
computer  modeling  has  shown  is  not  conducive  to  obtaining  initial 
distributions.  Therefore,  this  system  was  studied  again  with  the 
modified  experimental  apparatus. 

The  results  reported  are  from  a  CM-PP  experiment  in  which  a  D2  flow 
rate  of  5.2  moles/sec  and  a  F2  flow  rate  of  14  moles/sec  was  used. 
These  were  the  highest  flows  that  could  be  used  without  entering  the 
high  flow  region.  As  in  the  H  +  F2  system,  vibration-rotation 
populations  were  obtained  from  a  nonlinear  least  squares  fit  to  equation 
IV-2  assuming  a  Boltzmann  distribution.  The  rotational  temperature 
obtained  from  this  fit  was  400  K,  the  same  as  for  the  H  +  F2  reaction. 

The  DF (v)  vibrational  population-time  profiles  and  DF(v)  ratio  plots 
are  shown  in  Figures  V-ll  and  V-12.  In  order  to  avoid  congestion  in 
these  figures,  every  other  vibrational  level  is  shown.  The  product 
energy  distribution  is  compared  to  other  measurements  in  Table  V-4.  The 
CM-PP  results  are  in  agreement  with  the  Measured  Relaxation  results13 
peaking  at  v=10.  The  results  obtained  by  the  fast  flow  method59  peak  at 
v=9.  Usually  when  a  distribution  peaks  at  a  higher  viorational  level, 
it  is  indicative  of  the  true  nascent  distribution  free  of  relaxation 
processes  since  vibrational  energy  transfer  populates  lower  levels. 

A  comparison  of  the  results  of  the  D  +  F2  reaction  with  the  H  +  F? 
reaction  shows  a  higher  fraction  of  the  available  energy  being  channeled 
into  vibration  for  the  former  reaction.  This  is  in  contrast  to  the  F  + 
(H2/d2)  system  where  no  isotope  effect  is  observed.  Jonathon13  has 
suggested  that  the  H  +  F2  reaction  requires  more  translational  energy 
than  the  D  +  F2  reaction  to  surmount  its  slightly  greater  activation 
barrier.  Thus  the  hydrogen  atom  would  approach  closer  to  the  fluorine 
molecule  than  the  deuterium  atom  before  the  repulsive  energy  is 
released.  This  would  result  in  more  energy  being  channeled  into  recoil 
of  the  products  and  less  into  vibrational  motion.  Another  possible 
explanation  could  be  due  to  the  differences  in  the  masses  of  H  and  D. 

As  the  repulsive  energy  between  the  F2  is  released,  more  of  the 
available  energy  is  channeled  into  vibrational  motion  when  the  mass  of 
the  atom  is  heavier. 
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FFb 
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1 

— 

— 

0.00  +/-  0 

2 

0.06 

— 

0.00  +/-  0. 

3 

0.10 

0.01 

0.00  +/-  0 

4 

0.12 

0.03 

0.00  +/-  0 

5 

0.18 

0.07 

0.20  +/-  0 

6 

0.18 

0.14 

0.36  +/-  0 

7 

0.49 

0.22 

0.55  +/-  0 

8 

0.70 

0.48 

0.28  +/-  0 

9 

1.00 

1.00 

1.00 

10 

1.  19 

0.73 

1.15  +/-  0 

11 

0.75 

0.51 

0.66  +/-  0 

12 

0.29 

0.15 

— 

13 

<0.10 

0.04 

— 

14 

<0.08 

— 

— 

<f  > 

V 

0.62 

0.64 

0.62 

Measured 

Relaxation  from 

reference  13. 

^Fast  Flow 

from  reference 

59. 

'This  work. 
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F  +  HBr 

The  HF(v)  produced  from  this  reaction  can  populate  up  to  vibrational 
level  4  (42.4  kcal/mole55  above  the  zero  point  energy)  since  the 
available  energy  is  51  kcal/mole&0.  This  reaction  is  included  in  this 
work  since  there  is  some  disagreement  as  to  the  "correct"  nascent 
product  vibrational  energy  distribution. 

In  1979,  Polanyi's  group^  studied  the  F  +  HBr  reaction  using  the 
Arrested  Relaxation  technique  and  observed  a  dependence  on  reagent  flow. 
They  concluded  that  a  relaxed  distribution  describes  the  initially 
formed  products  in  the  F/HBr  system  and  proposed  that  at  high  flows 
preferential  deactivation  of  the  low  vibration  levels  cause  the 
distribution  to  become  inverted.  They  suggested  that  the  deactivation 
process  could  be  vibration-electronic  energy  transfer  by  Br  atoms. 

Subsequently,  Jonathon  did  an  arrested  relaxation  experiment  of 
the  F/HBr  system  in  which  he  added  excess  Br?.  Flourine  atoms  react 
with  Br2  three  times  faster  than  with  HBr  and  produce  ground  state  Br 
atoms  and  FBr.  2  He  observed  no  change  in  the  initial  distribution  with 
and  without  the  added  Br?,  indicating  that  deactivation  by  Br  atoms  was 
not  responsible  for  the  flow  dependence.  He  therefore  concluded  that 
the  initial  distribution  was  inverted  and  in  agreement  with  his  previous 
studies  of  this  reaction  in  which  no  flow  dependence  had  been  observed. 
Jonathon's  Arrested  Relaxation  technique  had  a  faster  pumping  system 
than  that  used  in  other  Arrested  Relaxation  techniques,  although  he  did 
do  some  experiments  in  which  the  pumping  speed  was  reduced  with  still  no 
change  in  the  observed  distributions. 

6  n 

Setser's  groupDU  obtained  results  using  an  Arrested  Relaxation 
technique  in  which  the  flow  dependence  was  observed  and  also  a  Fast  Flow 
technique  in  which  the  flow  dependence  was  absent.  The  distributions 
obtained  under  moderate  reagent  flow  conditions  from  the  Arrested 
Relaxation  experiments  were  inverted  and  in  agreement  with  the  results 
of  the  Fast  Flow  experiments.  He  favored  these  distributions  as  being 
the  nascent  product  distributions  and  suggested  that  wall  effects  were 
important  at  low  flows  in  the  Arrested  Relaxation  experiments.  Excited 
HF  molecules  could  reenter  the  observation  zone  after  partial  relaxation 
by  the  liquid  nitrogen  cooled  walls  or  the  HF  molecules  could  be  formed 
outside  the  observation  zone  and  then  reenter  after  collisional  or 
radiative  relaxation.  The  existence  of  some  type  of  wall  effect  was 
supported  by  the  fact  that  the  flow  dependence  was  not  observed  for  non¬ 
condensibles  such  as  H2,  CH4,  and  C2Hg. 

Sloan  recently  studied  the  F  +  HBr  reaction  using  a  low  pressure 
infrared  chemiluminescence  apparatus  in  which  three  concentric  quartz 
tubes  were  used  for  reagent  injection.  The  innermost  tube  could  be 
raised  or  lowered  thus  enabling  the  reagent  mixing  zone  to  be  moved  to 
regions  of  high  or  low  reagent  density  without  changing  reagent  flows. 
Thus  reaction  could  be  studied  over  a  range  of  conditions.  For  the 
experimental  configuration  in  which  deactivation  processes  were 
eliminated,  he  obtained  an  inverted  vibrational  product  distribution. 

He  also  showed  that  by  changing  reagent  flows  and  the  mixing  geometry 
that  vibrational  relaxation  by  the  walls  of  the  vessel  or  the  injector 
is  important  when  the  mean  free  path  in  the  reaction  vessel  is  unusually 
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large  as  it  is  for  fast  reactions  such  as  F  +  HBr.  He  pointed  out, 
however,  that  these  conclusions  were  obtained  by  operating  his  apparatus 
in  a  way  that  led  to  this  result;  not  the  normal  operating  condition. 

Table  V-5  compares  the  distributions  obtained  for  the  F  +  HBr 
reaction  from  experimental  methods  in  which  the  flow  dependence  was  not 
observed.  Included  in  this  table  are  the  results  obtained  from  the  CM- 
PP  and  CM-LP  experiments  of  the  present  work.  The  flows  used  in  the  CM- 
PP  experiment  were  CF^  =  2.6  mole/sec,  and  HBr  =  3.0  and  4.4  mole/sec. 
These  were  the  limiting  flows  possible  to  give  measurable  intensity  and 
still  not  to  be  in  the  high  flow  region.  The  presssures  used  in  the  CM- 
LP  experiments  were:  SFg  =  60  mtorr,  HBr=  10-50  mtorr,  and  Ar=  3.5  torr. 
Within  experimental  error,  the  distributions  obtained  from  both  the  CM- 
PP  and  CM-LP  experiments  showed  no  dependence  on  reagent  flow.  The 
rotational  temperatures  calculated  for  the  CM-PP  experiments  were  300  K. 
The  temperatures  of  the  distributions  obtained  from  the  CM-LP 
experiments  were  300  K  for  HF(v=4),  400  K  for  HF(v=3)  and  HF(v=2),  and 
500  K  for  HF(v=l). 

The  distributions  obtained  from  the  CM-PP  and  CM-LP  experiments  are 
in  closer  agreement  with  the  inverted  distributions  than  the  reported 
relaxed  distributions  obtained  in  experiments  where  a  flow  dependence  is 
observed.  The  CM-LP  distribution  is  slightly  broader  than  the  CM-PP 
distribution  lying  just  within  the  limits  of  experimental  error. 

Sloan's  group63  reports  the  most  inverted  distribution;  by  comparison 
all  of  the  other  distributions  show  some  degree  of  relaxation. 

In  the  present  CM-pp  and  CM-LP  experiments,  HF(v)  collisional 
relaxation  by  HBr,  HF(v*=0)  or  Ar  and  HF(v)  radiative  relaxation  can  be 
eliminated  as  causes  of  the  slightly  relaxed  distributions.  Computer 
simulations  have  shown  that  these  processes  do  not  affect  the 
determination  of  the  initial  distributions  under  the  conditions  used  in 
the  present  experiments.  Relaxation  by  the  walls  of  the  reaction  vessel 
in  the  CM-PP  experiments  is  unlikely  since  it  would  take  a  minimum  of 
about  0.5  msec  for  the  molecules  to  travel  to  the  wall  and  back  into  the 
observation  volume.  Since  the  extrapolation  is  performed  for  data  only 
out  to  0.5  msec,  the  distribution  should  not  be  affected. 

Table  V-6  is  a  comparison  of  the  reaction  and  deactivation  rate 
constants  (HBr  as  deactivator)  for  this  system.  The  reaction  rate 
constant  obtained  from  analysis  of  the  total  intensity  time  profile  is 
in  agreement  with  total  intensity  measurements  reported  by  Setser60  and 
Wurzburg  and  Houston64.  The  reaction  rate  constant  obtained  in  this 
work  from  the  v=4  time  profile  is  slightly  outside  the  error  limits  of 
the  value  reported  by  Smith  and  Wrigley65  who  analyzed  overtone  emission 
from  v=4.  Analysis  of  the  experimentally  observed  HF(v=3)  population¬ 
time  profiles  after  subtraction  of  the  relaxation  contribution  from 
HF(v=4)  gives  a  total  reaction  rate  constant  of  2.4  +/-  0.2  x  106  torr-1 
sec,-1  in  agreement  with  that  obtained  from  HF(v=4)  as  expected.  As  has 
been  shown  by  modeling  in  this  work,  total  intensity  profiles  give  low 
total  rate  constants  when  forced  to  fit  the  sum  of  two  exponentials, 
therefore,  the  result  from  the  fit  of  the  highest  vibrational  level  is 
preferred.  The  intercept  of  the  pseudo  first  order  reaction  rate 
constant  versus  HBr  pressure  was  zero. 
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Table  V-5 

Relative  Vibration  Distributions  for 


the 

Reaction  F 

+  HBr 

— ->  HF(v)  + 

Br 

V 

FFa 

FFb 

JLSl 

CM-gPe 

CM-LP® 

1 

0.26 

0.28 

0.35 

0.54 

0.45  +/-  0.22 

0  50 

+/-  0. 16 

2 

0.65 

0.67 

0.58 

0.68 

0.66  +/-  0. 10 

0.82 

+/-  0.10 

3 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4 

1.03 

1.11 

1.29 

1.05 

0.99  V-  0.18 

0.79 

+/-  0.06 

<v 

0.60 

0.63 

0.56 

0.60 

0.57 

<v 

0.13 

0.08 

0.11 

<v 

0.28 

0.30 

0.33 

aFrom  reference  60. 
bFrom  reference  51. 

cFrom  reference  63,  infrared  chemiluminescence  experiment. 
dFrom  reference  62,  arrested  relaxation. 
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Table  V-6 

Reaction  and  Deactivation  Rate  Constants 
for  the  Reaction  F  +  HBr  - >  HF(v)  +  Br 


Reaction  Rata  Conatanta 

( x  10®  torr'1  aac'1) 

*tot 

Mat,  had 

1.46  +/-  0.13 

f  froa  Photolraia  SR. 

HF( total  lntanaltr)a° 

1.5 

Taat  Flow6 

2.0  ♦/-  0.2 

F  froa  Photolrala  F, 

HF(v=4)  Tlaa  Raaolvtd  IR 

Chaal  luainaaconea0 

1.5  ♦/-  0.8 

d 

CM-LF  HF(  total  intanaitr) 

2.8  ♦/-  0.2 

CM-LP  8F(»a3)<1 

2.4  ♦/-  0.2 

CH-LP  BF(ys2)d 

Daactlvatlon  Rata  Conatanta 

ky(H8r)  s  10*  torr'1  aac’1 

v  =  1 

2  3 

4 

0.27  ♦/-  0.04* 

72  ♦/*  0 . 7C 

0.27  ♦/-  0.04* 

2.3  ♦/-  0.5d 

3.2  */-  0. 9d 

aFrom  reference  64.  eFrom  reference  71. 

^From  reference  60.  ^From  reference  72. 


cFrom  reference  65- 
^This  work. 


Rate  constants  for  deactivation  by  HBr  are  also  shown  in  Table  V-6. 
The  rate  constant  for  deactivation  of  HF(v=4)  by  HBr  obtained  in  this 
work  is  approximately  a  factor  of  2.3  lower  than  that  obtained  by  Smith 
and  Wrigley6  .  A  lower  deactivation  rate  constant  usually  means  that 
the  measurement  is  free  of  additional  deactivators.  Plots  of  pseudo 
first  order  deactivation  rate  constant  versus  HBr  pressure  gave  an 
intercept  of  3.3  +/-  0.9  x  103  sec""1.  As  contrasted  to  the  H2  system, 
the  relative  error  in  the  determination  of  the  slope  for  this  system  is 
not  significantly  affected  by  the  large  intercept  which  is  measured. 

The  relaxation  rate  constant  measured  for  HF(v=3)  deactivation  by 
HBr  is  2.3  +/-  0.5  x  105  torr-1  sec-1.  Measurements  of  this  rate 
constant  have  not  been  previously  reported.  Relaxation  rate  constants 
for  HF(v=2)  relaxation  by  HBr  could  not  be  extracted  from  the 
experimental  data  since  the  contribution  to  this  population  from  HF(v=3) 
relaxation  was  greater  than  70%.  When  the  total  intensity-time  profile 
was  fit  to  the  sum  of  two  exponentials,  a  deactivation  rate  constant  of 
3.9  +/-  0.3  x  104  torr  1  sec-1  was  obtained. 

F  +  CH2C12 

The  total  energy  available  to  this  system  is  44.7  kcal/mole66 
resulting  in  populated  vibrational  levels  up  to  v=4  (42.4  kcal/mole55 
above  the  zero  point  energy).  A  study  of  this  system  was  done  since  as 
with  the  F  +  HBr  reaction,  the  initial  distributions  obtained  from 
Arrested  Relaxation  experiments66'67  show  a  dependence  on  reagent  flow. 
In  addition,  deactivation  of  HF  by  CH2C12  has  not  been  previously 
reported. 

Polanyi's  group67  observed  a  flow  dependence  in  the  reactions  of  F 
atoms  with  the  substituted  methane  series  CH3CI,  CH2C12,  CHCI3.  As  in 
their  study  of  the  F/HBr  system,  they  attributed  the  results  at  high 
reagent  flows  to  be  due  to  preferential  deactivation  of  the  low 
vibrational  levels  and  reported  a  relaxed  distribution  as  being  the 
nascent  product  energy  distribution. 

A  collaborative66  effort  by  Setser's  group  in  Kansas  and  Heydtmann's 
group  in  Frankfurt  studied  the  same  methyl  chloride  series  as  Polanyi's 
group67  in  Toronto  and  a  comparison  between  the  results  obtained  in  the 
three  different  laboratories  was  made.  Quantitative  agreement  among  the 
three  laboratories  for  the  results  from  arrested  relaxation  methods,  in 
which  the  flow  dependence  was  observed,  was  poor.  For  the  F/CH2C12 
system,  only  distributions  favoring  the  v  =  l  level  were  observecT in 
Frankfurt.  Almost  all  of  the  arrested  relaxation  results  showed  a 
minimum  at  v=2  which  suggested  the  superposition  of  two  distributions; 
one  peaking  at  v  =  l  and  the  other  at  v  =  3.  When  F2  replaced  CF4  as  the 
source  of  F  atoms  in  the  Frankfurt  laboratory,  no  effect  on  the 
distribution  was  observed.  The  deviation  from  the  inverted  distribution 
of  Fast  Flow  measurments  was  found  to  be  greatest  for  CHCI3  which  has 
the  highest  number  of  chlorine  substituents  and  the  smallest  reaction 
rate  constant  of  the  series.  To  explain  the  arrested  relaxation  results, 
they  hypothesize  the  existence  of  two  reaction  channels.  One  channel 
produces  vibrational ly  excited  HF(v)  while  the  other  forms  a  metastable 
complex,  CH2C1o*F,  that  upon  dissociating  would  give  a  statistical 
distribution.  They  suggest  that  because  the  arrested  relaxation 


experiments  are  done  with  liquid  nitrogen  cooled  walls,  the 
translational  temperature  of  the  reagents  may  be  lower  than  300  K  so 
complex  formation  is  favored.  They  also  suggest  that  for  reactions 
which  have  small  rate  constants,  the  probability  that  the  reagent 
temperature  will  be  reduced  prior  to  reaction  is  greater.  In  addition, 
complex  formation  would  be  more  favorable  with  increasing  chlorine 
substitution. 

The  flows  used  in  the  CM-PP  experiments  were  CF/  =  3.7  mole/sec  and 
CH2CI0  =  1.9-3. 7  mole/sec.  The  pressures  used  in  tne  CM-LP  experiments 
were  SFg  =  60  mtorr,  CH2CI2  =  33  and  53  mtorr.  No  flow  dependence  was 
observed  in  either  experimental  method.  The  temperature  of  the 
rotational  distributions  from  the  CM-PP  experiments  was  300  K.  The 
temperatures  of  the  rotational  distributions  from  the  CM-LP  experiments 
were  500  K  for  the  first  vibrational  level,  400  K  for  the  second 
vibrational  level,  and  300  K  for  the  third  vibrational  level. 

Table  V-7  gives  the  nascent  product  energy  distributions.  The 
distributions  obtained  from  both  the  CM-PP  and  CM-LP  experiments  support 
the  inverted  distribution  obtained  from  the  Fast  Flow  experiments68  in 
which  no  flow  dependence  was  observed.  To  test  the  metastable  complex 
hypothesis,  CM-PP  and  CM-LP  experiments  could  be  done  in  which  the 
translational  energy  of  the  reagents  is  varied. 

The  reaction  rate  constants  for  the  CH2Cl2  system  are  given  in  Table 
V-8.  The  plot  of  the  pseudo  first  order  reaction  rate  constant  versus 
CH2C1 2  presssure  gave  an  intercept  of  zero.  The  reaction  rate  constant 
obtained  from  the  fit  of  total  intensity  is  57%  -  160%  larger  than  the 
previously  reported  reaction  rate  constants.  There  could  be  two  causes 
for  this  observation.  The  present  result  is  a  measure  of  F  atom 
depletion  so  that  if  some  other  species  present  is  depleting  F  atoms, 
this  would  result  in  an  apparent  faster  reaction  rate  constant. 
Alternatively,  F  could  be  depleted  by  complexing  with  CH2C12  or  reacting 
with  CHCI2  as  a  secondary  reaction. 

The  previously  reported  reaction  rate  constants  could  be  low.  In 
the  fast  flow  experiments68,  the  reaction  rate  constant  is  obtained  from 
the  measurement  of  the  formation  of  HF(v  >  0)  and  so  if  there  is 
appreciable  population  in  vibrational  level  0  then  the  reaction  rate 
constant  would  be  too  small.  The  other  two  reported  values  of  the 
reaction  rate  constant  measure  depletion  of  CH2CI2.  If  in  these 
experiments,  the  F  atoms  are  depleted  by  a  source  in  addition  to  CH2CI2 
then  these  reported  reaction  rate  constants  would  also  be  low. 

The  deactivation  rate  constant  obtained  from  the  analysis  of  the 
total  intensity  time  profile  was  5.2  +/-  0.5  x  104  torr-1  sec-1  and  the 
intercept  obtained  from  the  pseudo  first  order  relaxation  rate  constant 
versus  CH2CI2  plot  was  1.5  +/-  0.5  x  103  sec-1.  Unfortunately,  because 
of  the  small  signal/noise  ratio  of  these  experiments,  no  meaningful 
reaction  and  deactivation  rate  constants  from  the  analysis  of  the 
individual  HF(v)  levels  could  be  obtained.  The  results  from  analysis  of 
the  total  intensity  time  profiles  are  lower  limits  of  the  reaction  and 
deactivation  rate  constants  for  this  system  as  shown  by  the  computer 
simulations  and  the  previous  results  for  the  F  +  H2  system  and  the  F  + 
HBr  system. 


Table  V- 7 

Relative  Vibration  Distributions  for  the 
Reaction  F  +  CH2C12  — >  HF(v)  +  CHC12 


V 

_ee! 

-JM=g Pb 

gM-fc?b 

1 

0.17 

0.19  +/-  0.14 

0.44  +/"  0.42 

2 

0.31 

0.56  V  0.12 

0.48  +/-  0.26 

3 

1.00 

1.00 

1.00 

4 

0.08 

0.09  +/-  0.04 

<*„> 

0.66 

0.62 

0.56 

aFrom  reference  73. 
bThis  work. 
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Table  V-8 

Reaction  Rate  Constants  for  the 
Reaction  F  +  CHgClg  — >  HF(v)  +  CHC12 


k.  x  10^  t orr  *  s«c  * 

tot 

Method 

4.6 

Feat  Flow* 

4.9  +/-  0.9 

Diffusion  Cloudb 

3.0 

Mesa  Spec,  of  CH2C12c 

7.7  ♦/-  0.2 

CM-LP  (total  intensity) d 

aFrom  reference  73. 
bFrom  reference  74 
cFrom  reference  75 


VI.  CONCLUSIONS 


Previously  experiments  were  performed  to  determine  the  microscopic 
rate  constants  for  the  formation  of  DF(v)  from  the  reactions  F  +  D2  and 
D  +  F2  using  a  new  technique  called  Chemiluminescence  Mapping  (CM)  which 
involves  both  spectral  and  time  resolution.  The  results  of  these 
experiments  were  questioned.  Consequently  this  project  was  undertaken 
to  determine  the  validity  of  the  CM  method  and  to  confirm  the  early  CM 
results. 

The  work  reported  here  has  shown  that  energy  disposal  experiments 
in  which  a  microwave  discharge  is  used  for  the  formation  of  the  atomic 
reactant  is  not  contaminated  with  products  from  ion  molecule  reactions. 
Although  these  results  hold  specifically  for  CM-PP  systems  they  are 
valid  for  other  experiments  in  which  a  microwave  discharge  is  used  for 
the  formation  of  atomic  species.  Using  the  double  floating  probe 
technique  the  ion  content  for  the  CM-PP  configuration  was  found  to  be 
less  than  1.8  x  10“°  of  the  gas  flowing  through  the  microwave  cavity. 

Both  experiments  and  computer  simulations  reported  here  revealed  and 
evaluated  two  limiting  types  of  CM  experiments:  CM-PP  in  which  there  is 
a  repetitive  pressure  pulse  of  the  diatomic  species  (this  corresponds  to 
the  early  CM  system)  and  CM-LP  in  which  the  atomic  species  are  formed 
repetitively  in  a  constant  pressure  environment.  Both  types  of  CM 
experiments  are  capable  of  producing  nascent  vibrational  energy 
distributions.  CM-LP  experiments  can  also  provide  detailed  information 
on  the  rate  constants  for  reaction  and  relaxation.  The  advantages  of 
the  CM  technique  over  the  Arrested  Relaxation  technique  includes  a  low 
duty  cycle  reagent  pulse  which  reduces  average  molar  flow  and  pumping 
requirements,  the  absence  of  a  wavelength  dependent  baseline,  and  both 
spectral  and  time  resolution.  Arrested  relaxation  techniques  can 
provide  relative  reaction  rate  constants  buc  do  not  give  relaxation 
information.  Because  the  CM  technique  provides  both  spectral  and  tine 
resolution,  energy  transfer  information  can  be  obtained  as  a  function  of 
deactivator  and  vibrational  level.  Essentially  any  deactivator  may  be 
used,  unlike  fast  flow  experiments  in  which  only  quenchers  that  do  not 
react  with  the  atomic  reactant  can  be  studied.  Direct  excitation 
methods  can  be  used  to  study  any  quenchers  that  do  not  absorb  at  the 
same  wavelength  as  the  molecules  being  excited,  but  the  formation  of 
high  vibrational  levels  via  multiquantum  absorption  has  low  absorption 
coefficients  and  these  types  of  experiments  do  not  give  any  information 
on  reaction  rate  constants. 

Computer  simulations  showed  that  in  order  to  obtain  nascent 
vibrational  energy  distributions,  CM-PP  experiments  should  be  performed 
under  low  flow  conditions  (i.e.  F  >>  H 2 )  and  for  observation  times  short 
enough  to  resolve  the  fastest  relaxation  process.  The  reaction  vessel 
must  be  evacuated  between  pulse  cycles  to  reduce  complications  (i.e. 
secondary  reactions,  collisional  and/or  radiative  relaxation)  that  car. 
arise  from  residual  reactant  and/or  product  species.  Based  on  these 
results,  improvements  were  made  on  the  CM-PP  experiments  which  were 
previously  reported.  These  improvements  included: 

i)  replacement  of  the  6  channel  signal  averager  by  one  with  2048 
channels,  increasing  the  time  resolution  from  5  milliseconds  to  1 
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microsecond. 


ii)  replacement  of  the  diatomic  reactant  injector  with  a  piezoelectric 
valve,  reducing  the  pulse  length  from  10  to  0.5  milliseconds. 

iii)  reduction  of  observation  volume  from  1500  to  10  cc;  a  concomittant 
reduction  on  S/N  was  experienced. 

Analysis  of  the  CM-LP  experiments  is  more  direct  than  the  CM-PP 
experiments  because  of  the  pseudo  first  order  experimental  conditions. 
Simulations  were  used  to  show  that  total  reaction  and  microscopic 
relaxation  rate  constants  could  be  obtained  by  fitting  total  vibrational 
population-time  data  to  the  sum  of  two  exponentials.  For  vibrational 
levels  below  the  maximum  vibrational  level  populated,  it  is  necessary  to 
first  decouple  input  due  to  relaxation  of  the  upper  levels  from  input 
due  to  reaction.  The  simulations  revealed  that  fitting  total  emission- 
intensity  time  profiles  to  the  sum  of  two  exponentials  gave  reaction  and 
relaxation  rate  constants  that  were  too  small  by  20%  and  35% 
respectively. 

Experimental  results  were  presented  for  the  F  +  H2(D2),  H(D)  +  F2»  F 
+  HBr  and  F  +  CH2CI2  chemical  reactions.  A  comparison  of  the  nascent 
vibrational  distribution  and  reaction  rate  constant  (analysis  of  the 
total  intensity-time  profiles)  for  the  F  +  H2  reaction  obtained  from  the 
CM-PP  and  CM-LP  experiments  with  previous  values  for  this  system  are 
within  experimental  error  and  provide  a  successful  benchmark  for  the 
methods.  The  reaction  rate  constant  obtained  from  analysis  of  HF(v=3) 
is  1.4  +/-  0.3  x  106  torr  1  sec~  .  This  value  is  preferred  over  the 
smaller  value  obtained  from  total  intensity  measurements.  Microscopic 
relaxation  rate  constant  measurements  obtained  in  this  work  were  larger 
than  previously  determined  values.  The  reason  for  this  discrepancy  is 
believed  to  be  due  to  the  presence  of  H2O  in  the  CM-LP  experiments. 

A  comparison  of  previous  CM-PP  results  with  that  of  the 
present  work  for  the  F  +  H2/D2  systems  are  within  experimental 
error.  Computer  modeling  has  shown,  however,  that  collisional 
and/or  radiative  relaxation  affects  the  nascent  distributions  for 
the  experimental  conditions  of  the  previous  CM  experiments. 

These  effects  may  have  been  within  the  experimental  error  {>  10%) 
and  thus  not  apparent  in  the  F  +  H2/D2  experimental  results. 

The  nascent  distributions  of  the  present  CM-PP  experiments  for  the 
H/D  +  F2  reactions  are  within  experimental  error  of  other  results  for 
these  systems,  but  disagree  with  the  previous  CM  measurements.  It  is 
believed  that  these  differences  can  be  reconciled  by  the  double  maxima 
in  the  HF(v)  time  profiles  (indicating  high  flow  conditions),  the  large 
observation  times  in  the  earlier  experiments  and  an  absorption  'hole' 
which  appears  on  the  liquid  nitrogen  cooled  filter  of  the  detector 
dewar.  Collisional  and  radiative  relaxation  effects  on  the  nascent 
distribution  of  these  reactions  were  more  apparent  than  in  the  F  +  H2/D2 
systems  since  the  rate  of  reaction  is  much  slower  (factors  of  11  and  17 
for  the  D  and  H  reactions  respectively). 

The  initial  distributions  reported  here  for  the  F  +  HBr  and  F  + 
CH2CI2  reactions  support  the  previously  observed  inverted  distributions. 
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Mo  dependence  on  reagent  flow  was  observed  in  either  the  CM-PP  or  CM-LP 
experiments.  Total  reaction  rate  constants  and  deactivation  rate 
constants  for  KF(v=3,4)  relaxation  by  HBr  were  reported.  The  EF(v  =  3) 
deactivation  rate  constant  has  not  been  measured  prior  to  this  work. 

The  HF(v=4)  deactivation  rate  constant  is  smaller  than  the  previously 
reported  value,  indicating  that  the  present  measurement  is  free  of 
additional  relaxation  processes.  Deactivation  of  HF(v=3,4)  by  HBr  is 
approximately  10  and  20  times  more  efficient  than  deactivation  by  H2* 

Due  to  the  large  errors  in  the  F  +  CH2Cl2  system,  microscopic 
deactivation  rate  constants  could  not  be  obtained.  The  relaxation  rate 
constant  obtained  from  the  total  intensity  measurement  is  a  lower  limit 
of  HF  ( v )  relaxation  by  C^C^- 

Future  efforts  should  be  focused  on  improvements  in  S/K  and  the 
elimination  of  the  large  intercepts  in  pseudo-first  order  rate  constant 
versus  HR  pressure  plots  which  are  apparently  due  to  H20  contamination 
in  the  flow  lines,  vacuum  system  and  reaction  vessel.  The  analysis  of 
the  raw  experimental  data  can  be  improved.  In  the  present  work,  a  least 
squares  fit  of  the  spectrum  for  each  time  slice  is  made  in  order  to 
extract  vibrational-rotational  populations;  this  involves  inversion  of 
up  to  a  150  by  150  matrix  for  each  time  slice.  An  analysis  using  fast 
Fourier  Transform  filtering  and  a  constrained  least  squares  analysis 
using  the  single  valued  decomposition  method  is  presently  being 
developed.  It  is  anticipated  that  with  these  improvements  the  computer 
time  for  data  processing  will  decrease  by  an  order  of  magnitude. 
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